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I.  PRODUCTION 


This  report  decribes  both  theoretical  and  experimental  studies  into 
methods  for  producing  and  storing  metastable  atoms  and  molecules,  i.e.,  spe¬ 
cies  with  relatively  long  radiative  (usually  dipole  forbidden)  lifetimes.  The 
motivation  for  the  present  study  is  due,  in  part,  to  the  need  to  improve  the 
state  of  the  art  in  metastable  production  techniques.  Specifically,  there  is 
a  desire  to  increase  both  the  absolute  density  and  metastable  fractional  com¬ 
position  of  these  sources.  To  dav  ,  there  has  beer  little  effort  specifically 
directed  towards  these  particular  goals.  Any  method  which  accomplishes  one 
and/or  both  of  these  goals  will  provide  an  important  experimental  tool  that 
can  be  used  to  improve  the  data  base  of  the  fundamental  properties  of  meta¬ 
stable  atoms  and  molecules,  e.g.,  rate  constants,  collisional  and  excitation 
cross  sections,  and  lifetimes.  The  accurate  evaluation  of  these  quantities  is 
needed  to  understand  the  detailed  behavior  of  a  diversity  of  plasma,  dis¬ 
charge,  flame  and  other  equilibrium  and  nonequilibrium  situations  in  which 
processes  involving  metastable  atoms  and  molecules  play  an  important  role. 

Ihis  report  is  composed  of  several  separate  parts.  Section  II  is  a  lit¬ 
erature  search  and  compilation  of  radiative  lifetimes  of  some  representative 
metastable  atoms  and  molecules.  Dominant  radiative  decay  mechanisms  for  these 
states  are  also  reported  in  Section  II.  Section  III  is  *!evoted  to  a  brief  re¬ 
view  of  some  metastable  production  techniques,  and  an  experimental  study  of  a 
simple,  basic  research  source  that  can  be  used  to  produce  metastable  rare  gas 
atoms  and  metastable  nitrogen  molecules.  Section  IV  describes  novel  tech¬ 
niques  for  producing  a  high  flux  and/or  high  purity  source  of  metastable  spe¬ 
cies.  As  far  as  we  know,  none  of  these  techniques  have  been  experimentally 
verified  although  the  calculations  contained  herein  suggest  improvements  in 
flux  and  purity  of  (for  some  cases)  many  orders  of  magnitude.  Fina’ly,  some 
initial  thoughts  on  the  difficult  problem  of  metastable  storage  techniques  are 
described  in  Section  V. 
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II.  LITERATURE  SURVEY  OF  RADIATIVE  LIFETIMES 


A  literature  search  was  undertaken  to  compile  radiative  lifetimes  and, 
in  many  cases,  radiative  decay  mechanisms  (electric  dipole,  etc.)  for 
metastable  states  (i.e.,  states  with  radiative  lifetimes  >  1  us)  in  a  number 
of  atomic  and  molecular  systems.  The  candidate  species  were  chosen  solely  on 
the  basis  of  their  being  representative  of  systems  that  are  extensively 
studied  in  atomic  and  molecular  physics.  The  results  of  this  phase  of  the 
program  are  reported  in  Table  1.  The  first  column  in  the  table  identifies  the 
species  of  interest,  the  second  column  designates  the  particular  metastable 
level,  the  third  column  designates  the  energy  of  that  level,  the  fourth  column 
indicates  the  range  of  measured  or  calculated  lifetimes  for  that  level,  the 
fifth  column  indicates  the  principal  radiative  decay  mechanism  (not 
necessarily  to  the  ground  state)  and  the  last  column  is  the  literature 
citation. 

No  attempt  has  been  made  here  at  either  a  complete  or  critical  review  of 
all  lifetime  data.  Only  a  limited  number  of  the  known  metastable  states  for 
each  species  was  included  although  we  obtained  as  much  data  as  possible  for 
each  candidate  state.  In  some  cases  we  have  taken  the  liberty  of  deleting 
lifetime  values  when  it  is  apparent  that  the  literature  suggests  this  value  is 
in  error.  In  a  number  of  other  cases,  unresolved  controversies  still  exist. 

In  those  cases,  we  have  listed  all  the  reported  lifetime  values  in  orier  to 
give  an  estimate  of  the  precision  of  each  lifetime  value. 


TABLE  1  (Continued) 


4->  <0 


V) 


lu 

CO 

<  I 
h-  I 
CO  * 

?! 

LU 

£ 


<U  ■ 
4J 


o 


II 

+  i 


w 

CO 

< 


D1 


U1 


65 

**• 

<0 


+ 

W 

CM 

< 


65 

CM 

CO 


tS 

CO 

<o 


CL 

•*“  o 

HP  111 

O  O  E 
aw-r 
E'rP 
O  ,0) 
U<—  0- 

OL-  r— 

"53 

OCMPJ  O 
CL.  f0  II 
CO  C  > 


•r-  <U  *U 

E  T» 

0) 

r-  4J  4J 

10  O  <0 

c  c  +■> 
o  u> 

•r-  -P 
4J  3  r- 
«JQ  ID 
■M  C 

0-0 
S-  >>T- 
f—  4J 

11  11  K) 
-X  4J 

s:  -r-  O 

DC  r—  L- 


00 

UJ 

o 

LU 

CL 

00 


10  X)  0*0  QJ  4- 


CM 

o 


CM 


o 

z 


o 

C_J 


00 

LU 

H 

O 

z 


TABLE  1. 


1.  Crake,  G.W.F.,  Phvs.  Rev.  A3,  908  (1971). 

2.  Rondel,  R.O.  and  Stebbings,  R.F.,  Case  Studies  in  At .  Coll.  Phys.  ?  ed. 
E.W.  McDaniel  and  M.R.C.  McDowell  (North-Hoi land  Publishing' CoV, 
Amsterdam,  1972),  p.  550. 

3.  Small-Warren,  N.E.  and  Chiu,  L.Y.C.,  Phys.  Rev.  All,  1777  (1975). 

<1.  van  Dyck,  Or.,  R.S.,  Johnson,  C.E.  and  Shugart,  H.A.,  Phys.  Rev.  A5,  991 
(1972). 

5.  Garstang,  R.H. ,  J.  Opt.  Soc.  Amer.  52,  845  (1962). 

6.  Furcinitti,  P.S.,  Wright,  0.0.  and  Balling,  L.C.,  Phys.  Rev.  A12,  1123 
(1975). 

7.  Garstang,  R.H. ,  Astrophys.  0.  148,  579  (1967). 

8.  Hyman,  H.A.  (unpublished). 

9.  Garstang,  R.H.,  J.  Res.  NBS  68A,  61  (1964). 

10.  Shemansky,  D.,  J.  Chem.  Phys.  51.,  689  (1969);  Zipf,  E.C.,  Can.  0.  Chem. 

47,  1863  (1969);  Shemansky,  D.,  and  Carleton,  N.P.,  0.  Chem.  Phys.  51, 
b82  (1969).  “ 

11.  Tilford,  S.G.,  Wilkinson,  P.G.  and  Vandc-rslice,  J.T.,  Astrophys.  0.  141, 
427  (1965). 

12.  Shemansky,  O.E. ,  0.  Chem.  Phys.  51,  5487  (1969). 

13.  Lichten,  W.,  0.  Chem.  Phys.  26,  306  (1957). 

14.  Olmstead  III,  0.,  Newton,  A.S.  and  Street,  Or.,  K.,  0.  Chem.  Phys.  42, 

2321  (1965).  “ 

15.  Borst,  W.L.  ana  Zipf,  E.C.,  Phys.  Rev.  A3,  979  (1971). 

16.  Kaslow,  D.E.,  Ph.O.  Thesis,  Univ.  of  Michigan,  Ann  Arbor,  Mich.  (1971). 

17.  Holland,  R.F.,  J.  Chem.  Phys.  51,  3940  (1969). 

18.  Freund,  R.S.,  0.  Chem.  Phys.  50,  3734  (1969). 


6 


19.  Til  ford,  S.G.  and  Benesch,  W.M.,  0.  Chem.  Phys.  64,  3371  (1976). 

20.  Badg|r,  R.H. ,  Wright,  A.C.,  Whitlock,  R.F.,  0.  Chem.  Phys.  43,  4345 

21.  Nichol Is,  R.W.,  Can.  0.  Chem.  47,  1847  (1969). 

21a.  Vallance-Jones,  A.  and  Gattinger,  R.L.,  Planet.  Space  Sci.  11,  961 
(1963). 

22.  Wallace,  L.  and  Hunten,  O.M.,  0.  Geophys.  Res.  73,  4813  (1968). 

23.  Allen,  C.W.,  Astrophys.  J.  85,  156  (1937). 

24.  van  de  Hulst,  H.C.,  Ann.  Astrophys.  8,  12  (1945). 

25.  Childs,  W.H.J.  and  Mecke,  R.,  Z.  Physik  68,  344  (1931). 

26.  de  Jager,  C.,  Bull.  Astron.  Inst.  Neth.  13,  9  (1956). 

27.  Jarmain,  W.R.  and  Nicholls,  R.W.,  Proc.  Phys.  Soc.  90,  545  (1967). 

28.  James,  H.M.  and  Coolidge  A.S.,  Astrophys.  J.  87,  438  (1938). 

29.  Nicholls,  R.W.,  Ann.  Geophys.  20,  144  (1964). 

29a.  Bubert,  H.  and  Froben,  F.W.,  Chem.  Phys.  Lett.  8,  242  (1971). 

29b.  Oeunehomme,  M.,  J.  Chem.  Phys.  45,  4433  (1966). 

29c.  Jeunehomme,  M.  and  Ouncan,  A.B.F.,  0.  Chem.  Phys.  41^,  1692  (1964). 

29d.  Fink,  E.H.  and  Welge,  K.H.,  Z.  Na turf orsch  23a,  358  (1968). 

30.  Johnson,  C.E.  and  van  Dyck.  Jr,.F  R.S.,  J.  Chem.  Phys.  56,  1506  (1972). 

31.  James,  T.C.,  J.  Chem.  Phys.  55,  4118  (1971). 


7 


III.  DEVELOPMENT  OF  A  BASIC  RESEARCH  SOURCE  OF 
METASTABLE  ATOMS  AND  MOLECULES 

A.  INTRODUCTION 

The  objective  of  this  experimental  phase  of  the  program  is  the  develop¬ 
ment  of  techniques  for  producing  sources  of  metastable  atomic  and  molecular 

8  -3 

beams  that  can  generate  metastable  species  at  densities  j>  M  cm .  In 
particular,  there  is  an  interest  in  sources  for  rare  gas  electronic  metastable 
atoms  and  metastable  vibrational  or  low  lying  electronic  states  in  molecules. 
It  is  also  desired  that  the  fraction  of  metastable  species  relative  to  the  to¬ 
tal  beam  content  be  as  large  as  possible.  The  relatively  large  absolute  den¬ 
sity  and  high  fractional  purity  are  desirable  characteristics  for  enhancing 
the  signaT-to-noise  ratio  in  various  types  of  metastable  beam  scattering  ex¬ 
periments.  These  experiments  are  usually  difficult  even  when  both  collision 
partners  are  in  the  ground  state.  Low  beam  densities  and  concomitant  low  sig¬ 
nal  levels  are  required  for  producing  the  well-collimated  beams  needed  for 
making  angular  measurements.  The  level  of  experimental  difficulty  is  then 
further  increased  when  information  on  excited  state-ground  state  or  excited 
state-excited  state  collision  cross  sections  is  desired.  Nevertheless,  this 
information  is  of  fundamental  importance  for  understanding  the  properties  of 
excited  states  of  atoms  and  molecules,  and  the  development  of  methods  to  in¬ 
crease  the  data  base  in  this  area  should  be  pursued. 

Ouring  the  initial  phase  of  the  program,  it  was  decided  that  a  single 
candidate  technique  that  could  generate  both  rare  gas  atomic  and  molecular 
metastable  states  would  be  sought  in  order  to  satisfy  both  the  technical  re¬ 
quirements  and  the  scope  of  the  program.  This  strategy  has,  however,  pre¬ 
vented  us  from  investigating  some  interesting  but  untested  concepts  which 
might  greatly  improve  the  state  of  the  art  in  metastable  sources.  Of  special 
note  is  the  optical  pumping  method  for  the  production  of  low  lying  metastable 
levels  in  the  heavier  rare  gas  atoms  (see  Section  IV).  This  technique  is 
believed  capable  of  producing  mecastable  fractional  excitation  >  10  percent 
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but  was  not  considered  as  a  candidate  method  because  of  its  complexity  and  in¬ 
ability  to  produce  molecular  metastables  in  its  present  embodiment.  The 
choice  of  supersonic  nozzle  sources  was  also  discounted  due  to  the  limited 
pumping  speed  of  the  available  vacuum  system. 

B.  METAS7ABLE  SOURCES  REVIEW 

We  briefly  review  a  few  of  the  metastable  production  concepts  that  were 
considered  for  this  phase  of  the  program.  Section  III-B-5  reviews  the  method 
of  choice,  namely,  a  low  pressure,  high  current  discharge  source.  A  descrip¬ 
tion  of  the  experimental  arrangement  of  the  discharge  source,  associated  diag¬ 
nostics  and  results  follows  in  Section  III-C  and  III-D. 

1.  Atomic  Metastable  Production  via  Ion  Impact  on  Surfaces 

Experiments  performed  several  years  ago'  '  demonstrated  that  when  en¬ 
ergetic  helium  positive  ions  collide  with  a  metal  surface,  metastable  states 
of  the  neutral  parent  species  are  produced.  Subsequently,  this  technique  was 
employed  for  producing  beams  of  metastable  inert  gas  atoms  for  studies  of  sec¬ 
ondary  electron  emission  produced  by  metastable  atom  impact  on  solid  metal 
targets.  Ihe  more  recent^33)  of  these  publications  appeared  in  1950,  the 
original  observations  were  published  in  1929. Little  or  no  quantitative 
information  describing  the  metastable  beam  densities  and  fractional  population 
ratios  was  available  from  these  publications  since  the  main  thrust  of  the  work 
dealt  with  the  energy  distribution  of  the  secondary  electron  ejected  from  the 
metal  surface  and  hence  this  information  was  not  critical  to  the  needs  of  the 
experiments. 

Hagstrum^3^  has  reported  the  reflection  coefficients  (or  probability) 
for  producing  metastable  atoms  via  rare  gas  ion  impact  on  surfaces  as  a  func¬ 
tion  of  incident  Ion  energy.  A  total  of  five  surfaces  were  investigated: 
tungsten,  molybdenum,  silicon  (100),  hafnium  and  germanium  (111).  The  results 

of  this  work  indicate  that  :  (1)  the  reflection  coefficient  increases  rapidly 

with  ion  energy;  (2)  the  magnitude  of  the  reflection  coefficient  as  a  function 
of  velocity  is  independent  of  the  mass  or  nature  of  the  ion  (He  ,  Ar  , 

311  Oliphant,  M.L.E.,  Proc.  Roy.  Soc.  A124,  288  (1929). 

33.  Green,  0.  Phys.  Rev.,  876  (1950). 

34.  Hagstrum,  H.O.,  Phys.  Rev.  123,  758  (1961). 
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Ne  );  and  (3)  the  magnitude  of  the  reflection  coefficient  depends  upon  the 

cleanliness  of  the  surface.  For  atomically  clean  tungsten  and  molybdenum,  the 

-3 

reflection  coefficient  increases  from  about  10  at  100  eV  to  ~  0.04  at 

3 

10  eV.  For  atomically  clean  silicon,  the  values  are  about  two  to  ten  times 
higher.  Much  larger  values  have  been  observed  for  contaminated  germanium  end 
hafnium  surfaces  for  which  the  reflection  coefficient  approaches  0.1  at 
1000  eV. 

In  order  to  take  advantage  of  the  surface  conversion  of  ions  to  meta¬ 
stables,  a  glancing  incident  collision  is  required  in  order  to  remove  the 
metastable  from  the  proximity  of  the  surface  and  thus  prevent  re-ionization  of 
the  metastable.  This  application  was  reported  in  a  post-deadline  paper  at  th^ 
Paris  ICPEAC  Conference^)  which  described  an  experimental  measurement  of 
electron  impact  excitiation  cross  sections  on  metastable  helium  atoms  produced 
by  neutral izaion  of  positive  ions  on  a  graphite  surface.  Ion  optics  were  used 
to  focus  the  ions  from  a  Penning  ion  source  through  an  orifice  drilled  in  the 
graphite  plate.  Few  other  experimental  details  were  available  although  meta- 
stable  helium  densities  of  10  cm  were  quoted  and  cross-section  mea¬ 
surements  were  presented  establishing  the  utility  of  this  technique.  Although 
our  follow-up  inquiries  to  the  Russian  authors  requesting  technical  details 
has  met  with  no  response,  we  believe  that  a  quantitative  reinvestigation  of 
this  technique  is  required  to  validate  the  quoted  metastable  production 
especially  at  low  translation  ion  energy. 

2.  Charge  Exchange  Method 

The  basic  idea  is  to  utilize  the  fact  that  the  charge-transfer  cross 
section  becomes  very  large  (>  50  A^)  for  a  near  resonant  process.  If  one 
can  therefore  identify  charge-transfer  partners  in  which  capture  into  an  ex¬ 
cited,  metastable  state  is  near- resonant,  one  can  hope  to  selectively  produce 
a  large  density  of  metastable  atoms.  The  physical  principle  behind  the  method 
is  illustrated  by  the  energy-level  diagram  in  Figure  1.  This  figure  shows  the 
case  where  the  energy  defect  (aE)  relative  to  excited-state  A*  is  much  smaller 
than  that  relative  to  the  ground  state  A.  Under  such  circumstances,  one  would 
expect  the  reaction 

Kaxhaev,  A.O.,  Gostev,  V.A. ,  Saits'ev,  Yu.V.,  International  Conference  on 
the  Physics  of  Electronic  and  Atomic  Collisions,  Paris,  (1977). 
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Figure  1 


A+  +  B  -»  A*+  B+-  AE 


G5707 

Energy  Level  Diagram  for  Charge  Exchange  into  Excited  States 
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+  + 

A  +  B»a*  +  B  -  aE 

to  be  the  dominant  charge-transfer  process.  The  situation  shown  schematically 
in  Figure  1  is  typical  of  rare  gas  ion-alkali  atom  pairs.  Total  cross  sec¬ 
tions  for  He*  and  Ar+  ions  charge-transferring  to  K,  Rb  and  Cs  atoms  have 
been  measured  by  Peterson  and  torents.^3^  They  found  cross  sections 
>  50  A2,  and  although  they  did  not  attempt  to  directly  detect  excited  spe¬ 
cies,  they  concluded  that  such  large  total  cross  sections  could  only  be 
compatible  with  electron  capture  into  excited  states  of  He  and  Ar.  Their  con¬ 
clusions  were  subsequently  supported  by  the  theoretical  results  of  Ice  and 
Olson, who  noted  in  addition  that  the  combinations  Kr+  +  Rb,  and 

Xe+  +  Cs  should  favor  production  of  metastable  states  of  Kr  and  Xe.  Fin- 

1381 

ally,  Neynaber  and  Magnuson'  '  directly  measured  the  composition  of  He,  Ne, 
and  Ar  beams  resulting  from  the  rare  gas  ion-alkali  atom  charge  transfer  proe¬ 
ms,  and  found  that  >  50  percent  of  the  rare  gas  atoms  were  in  excited  meta¬ 
stable  states. 

A  conceptual  design  for  a  charge  transfer  source  is  shown  in  Figure  2. 
lhe  duoplasmatron  is  a  dc  arc-discharge  source,  yielding  a  high-intensity 
(-  1  mA),  space-charge-limited  beam  of  positive  ions  in  the  kilovolt  energy 
range.  The  electrostatic  Einzel  lens  is  a  focusing  device,  which  can  produce 
a  focused  beam  of  ions  -  0.1  cm  in  diameter  over  a  large  range  of  distances, 
lhe  beam  would  be  focused  in  the  interaction  region  (where,  for  example, 
electron-metastable  excitation  or  ionization  experiments  could  be  carried 
out),  which  would  be  separated  from  the  ion  source  chamber  by  a  small  aperture 
and  by  differential  pumping.  Steering  plates  would  permit  precise  alignment 
of  the  beam.  Prior  to  entering  the  interaction  region,  the  ion  beam  would 
transverse  a  well-defined  perpendicular  source  of  appropriate  alkali  vapor, 
effusing  from  a  narrow  slit  in  the  heated  oven.  The  perpendicular  geometry 
minimizes  the  amount  of  alkali  vapor  component  in  the  main  beam. 

TTbl  Peterson,  J.R.,  and  Lorents,  D.C.,  Phys.  Rev.  182,  152  (1969). 

37.  Ice,  G.E.,  and  Olson,  R.E.,  Phys.  Rev.  All,  111  (1975). 

38.  Neynaber,  R.H.,  and  Magnuson,  G.D.,  J.  Chem.  Phys.  65,  5239  (1976). 
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Figure  2  Conceptual  Design  for  a  Charge  Exchange  Metastable  Source 
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The  charge  exchange  technique  has  several  advantages  compared  with  other 
conventional  methods.  Remote  positioning  of  the  source  is  possible,  since  the 
beam  can  be  focused  as  required.  Since  the  cross  sections  are  very  large, 
conditions  in  the  oven  can  be  conveniently  arranged  to  guarantee  essentially 
complete  neutrality  of  the  emerging  beam.  Furthermore,  the  beam  will  contain 
~  50  percent  metastable  atoms,  whereas  for  conventional  methods  the  ratio  of 
metastables  to  ground-state  atoms  is  ~  10”3  -  10”^.  Finally,  the  absolute 
number  density  of  metastables  attainable  with  such  a  source  is  estimated  to  be 
~  10***  cm”'*.  The  charge-transfer  technique  for  producing  metstables 
should  certainly  be  useful  for  all  of  the  rare  gases  and  perhaps  also  for 
other  species  for  which  the  appropriate  charge-exchange  partners  can  be  found. 

3.  Optical  Pumping  of  Metastable  Calcium  Atoms 

The  production  of  the  lowest  lying  metastable  triplet  level  in  calcium, 
4s4p  at  1.88  eV,  has  been  accomplished  by  direct  optical  excitation  of 
the  intercombination  line  (4s^  *S  >4s4p  '^P1 )  at  6572.8  A  using  a 
tuned  dye  laser  to  excite  calcium  in  an  oven  at  ~  1100  K.v  '  The  oscil¬ 
lator  strength  of  this  transition  is  large  enough  that  a  significant  fraction 
of  the  ground  state  population  can  be  pumped  into  the  metastable  level  during 

p 

the  intense  dye  laser  pulse  (fluence  =  10  mJ/cm  ;  pulse  width  =  30  ns).  The 

calcium  density  was  controlled  by  varying  the  oven  temperature.  At  low 

optical  depths  corresponding  to  calcium  densities  of  the  order  of 

3  x  10* 3  cm”3  to  b  x  lu*4  cm3,  the  authors  report  about  50  percent  of 

the  atoms  excited  in  a  cylinder  volume  (1  cm  by  10  cm  long).  At  high 

optical  depths,  ~  10  cm”’  calcium  atoms,  the  absolute  density  of  excited 

atoms  can  be  obtained  from  a  measurement  of  the  laser  energy  since  the  number 

of  photons  in  the  laser  beam  is  much  less  than  the  number  of  calcium  atoms  in 

the  line  of  sight.  From  the  measured  laser  energy  (~  7  iml),  a  measurement  of 

the  line  width  of  the  intercombination  line  (~  0.15  A)  and  a  measurement  of 

the  laser  line  width  (~  0.5  A),  the  authors  obtain  an  absolute  metastable 
15  -3 

density  of  ~  1U  cm  under  these  condtions. 


3gT  Mcllrath,  T.J.,  and  Carlsten,  J.C.,  0.  Phys.  B6,  697  (1973). 
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4 .  Production  of  02  (*Ag)  via  Chemical  and  Photolytic  Methods 
Solution  phase  reactions  of  the  type: 


(?) 


have  long  been  known  to  produce  the  metastable  oxygen  molecule  O^Ag).  The 

mechanism  was  originally  elucidated  by  Kasha  and  co-workers. Recently, 

this  reaction  has  been  used  to  provide  a  chemical  pumping  mechanism  for  the 

0o/ 1  chemical  iodine  laser.  Measurements  of  the  absolute  0o(*a  1  and 
c  3  _  c  9 

02(  Eg")  ground  state  density  generated  in  the  reaction  of  alkaline 

Ho09  and  OC  “  were  performed  in  order  to  determine  whether  sufficient 

c  c  Mil 

excited  oxygen  is  produced  to  produce  iodine  atom  lasing.'  *  All  species 

concentrations  were  obtained  by  an  electron  paramagnetic  resonance  technique 
using  ground  state  02  as  a  calibrant  and  calculated  matrix  elements  to  de¬ 
termine  the  relative  signal  strengths  of  the  other  species.  At  the  highest 
Cu  flow  rate  used,  13  cm  /sec,  about  10  percent  of  the  oxygen  is  produced 

.  c  i 

in  the  A^state.  The  absolute  metastable  density  was  reported  to  be 
~  10  3  cm  .  The  90  percent  of  the  02  produced  in  the  ground  state  was 
attributed  to  the  base  catalyzed  decomposition  of  HgO,,. 

Another  unique  02(1Ag)  production  technique  utilizes  photolysis  of 
0^  below  about  300  nm  via  the  process: 

03  +  hv  (x  <  310  nm)  >  0(10)  +  O^Ag)  (3) 


which  has  an  energetic  threshold  near  310  nm.  The  relative  quantum  yield  for 
0(10)  production  has  been  meaiured  to  be  1.0  *  0.1  in  the  region  250-290  nm. ^42) 
The  absolute  quantum  efficiency  for  0(*D)  production  has  been  measured 

3TT  Kahn.  A.U. ,  and  Kasha,  M.,  0.  Amer,  Chem.  Soc.  92,  3293  (1970). 

41.  McDermott,  W.E.,  Lotz,  R.W.,  Delong,  M.L.,  and  Thomas,  D.M.,  "Chemical 
Generation  of  Og^Ag)". 

42.  Fairchild,  P.W.,  and  Lee,  E.K.C.,  Chem.  Phys.  Lett.  60,  36  (1978). 
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(43) 

to  be  1.0  at  253.7  nm  for  0^  immersed  in  liquid  Ar.  '  We  are  not  aware 

of  quant  r,  y'eld  measurements  for  09(*a  )  production.  However,  we  would 

1  C  3 

expect  a  significant  0^(  Ag)  production  via  KrF  (,\  ■  249  nm)  laser  pho¬ 
tolysis  of  0,  assuming  that  the  quantum  y  elds  of  0{*D)  and  09{V_) 

J  -17  z  9 

are  similar.  Ihe  03  absorption  coefficient  is  about  10  cm  at  Us 

peak  near  249  nm.  This  value  implies  that  a  laser  fluence  of  the  order  of 
0.1  J/rm  would  be  required  to  insure  a  high  photolysis  probability  per  in¬ 
cident  photon.  A  total  of  ~  10^  0g ( 1  )  could  be  generated  using  a 
100  mj  KrF  laser  pulse.  Both  the  fluence  and  laser  power  are  readily  attain¬ 
able  using  commercial  KrF  lasers. 

5.  Low  Pressure,  High  Current  Discharge 

After  reviewing  a  number  of  possible  alternatives,  it  was  decided  that  a 
low  pressure,  high  current  discharge  of  the  type  originally  described  by  sev¬ 
eral  groups^ was  the  method  of  choice  for  the  experimental  program  de¬ 
scribed  below.  These  reports  indicate  that  metastable  beams  of  alkaline  earth 
atoms,  e.g.,  Ca(3PQ  ^  ?),  Ca(*Dg),  Mg(3PQ  j  2)  and  Ba(30^  2  3)»  can  be  pro¬ 
duced  with  fractional  excitation  levels  approaching  30  percent  using  high  dis- 

o 

charge  current  fluxes  of  the  order  of  1  A/cnT  and  low  atomic  vapor  pressures, 

•aiS 

a  2  x  10  torr,  where  the  source  is  operating  under  effusive  conditions. 

The  low  voltage  (~  100  V),  high  current  discharge  was  obtained  by  accelerating 
thermionically  generated  electrons  in  front  of  the  exit  orifice  of  e  resis- 
tively  heated  oven  containing  tne  appropriate  aixaline  earth  atomic  vapor. 

Although  a  detailed  model  of  the  source  operation  is  not  available, 
there  are  several  qualitative  aspects  which  suggest  that  this  source  may  also 
be  applied  to  the  production  of  metastable  species  of  interest  here.  First, 
the  high  current  flux  is  necessary  in  order  to  obtain  high  metastable  produc¬ 
tion  rates  via  direct  electron  impact  excitation  of  ground  state  species: 


es  +  M  >  M*  +  e^  (4) 

DeMore,  W.B., and  Raper,  O.F.,  J.  Chem.  Phys.  44,  1780  (1966). 

44.  Brinkmann,  U.,  Goschler,  J.,  Steudel,  A.,  and  Walther,  H.,  Zeits.  fur 
Physik  228,  427  (1969);  Brinkmann,  U.,  Steudel,  and  H.  Walther,  Z.  Angew, 
Phys.  2fT223  (1967). 

45.  Isliii,  S.,  and  Qhlendorf,  W.,  Rev.  Sci.  inst.  43,  1632  (1972). 
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where  es  denotes  a  discharge  electron;  M  and  M  denote,  respectively,  the 
ground  and  metastable  levels;  and  the  energy  threshold  of  this  process  is 
given  by  E(e  ')  -  £(e,).  The  M  production  rate  is  linear  in  discharge 
electron  density  up  until  M  loss  processes  become  important  at  high  M 

if 

density  at  which  point  the  M  density  begins  to  level  off.  Second,  a  low 
atomic  beam  density  is  necessary  to  minimize  heavy  particle-metastable  loss 
processes.  Also,  if  the  atomic  beam  density  is  low,  the  electron  energy  loss 
can  be  made  small  enough  so  that  the  electron  average  energy  is  directly  de¬ 
termined  by  the  applied  field  strength. 

In  order  to  obtain  a  high  electron  current,  some  fraction  of  the  atomic 
beam  must  be  ionized  to  neutralize  the  negative  space  charge  which  tends  to 
limit  the  amount  of  electron  current  that  can  be  drawn  from  an  electron 
source.  For  plane  parallel  geometry,  this  limit  is  given  by  the  Chi  1 d- 
Langmuir  law:^6^ 


3/2 

J_.v  *  2.33  x  10“6  ^-5-  amp/cm^  (5) 

IMuX  Cm 


where  0  is  the  maximum  electron  current  flux  which  can  be  drawn  across  a 
region  of  space  bounded  by  two  electrodes  having  a  voltage  difference  V  be¬ 
tween  them  and  spaced  d  centimeters  apart.  This  is  about  the  limiting  current 
that  has  been  observed  in  space-charge  limited  electron  bombardment  metastable 

sources. For  typical  cases:  d  «  1  cm,  V  =  30  volts,  0m.v  «  385  yamps/cm^ 

M  7  3  HiaX 

and  the  maximum  electron  density,  n„,  -  10'  cm  assuming  an  average 
electron  energy  equals  30  eV.  This  value  of  n^j  is  smaller  than  that 
achievable  with  partial  space  charge  neutralization. 
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The  requirements  of  high  electron  current  (density)  and  low  neutral  beam 
density  cannot  be  achieved  simultaneously  under  all  conditions.  On  the  one 
hand,  the  neutral  beam  density  must  be  high  enough  to  maintain  a  sufficient 
ionization  rate  needed  for  charge  neutralization.  On  the  other  hand,  the  neu- 
t  .1  beam  density  should  be  low  enough  to  minimize  heavy  particle  metastable 
qL  nching.  In  addition  to  beam  density,  the  voltage  across  the  disch,  rge  is 
another  parameter  which  must  be  optimized.  Metastable  excitation  cross  sec¬ 
tions  typically  peak  at  much  lower  energies  than  ionization  cross  sections 
implying  that  the  ion  production  rate  may  not  maximize  at  the  same  applied 
voltage  for  which  the  metastable  production  rate  maximizes.  In  the  studies 
described  below  we  have  investigated  the  performance  of  the  high  current  dis¬ 
charge  source  as  a  function  of  discharge  current,  voltage  and  source  pressure 
in  order  to  optimize  the  metastable  output  for  argon,  krypton  and  N2  beams. 

C.  APPARATUS 

The  experiment  was  housed  within  a  differentially  pumped  stainless  steel 
vacuum  apparatus  shown  in  Figure  3.  The  metastable  source  and  detection  sys¬ 
tem  were  each  located  in  separate  sections  of  the  apparatus  and  communicate 
with  each  other  via  a  0.038  cm  diameter  orifice  (see  Figure  4).  Each  section 
is  pumped  by  a  6  in.  oil  diffusion  pump  and  can  attain  ~  10-®  torr  ultimate 
vacuum.  The  quoted  net  pumping  speed  of  the  diffusion  pump,  gate  val  /e  and 
liquid  nitrogen  trap  combination  stack  is  about  800  1/sec  (air). 

The  discharge  source  was  located  in  the  source  chamber  (Figure  4)  and 
was  comprised  of  a  pair  of  stainless  steel  bias  plates,  a  quartz  nozzle  gas 
source  and  a  heated  thoria  coated  iridium  filament.  A  picture  of  t  e  dis¬ 
charge  source  is  shown  in  Figure  5  and  •  schematic  of  the  discharge  source  and 
detector  is  shown  in  Figure  6.  The  gas  nozzle  source  was  made  by  h  ating  6  mm 
quartz  tubing  to  white  heat  and  then  drawing  the  tubing  apart  until  thin 
(0.018  cm  diameter)  capillary  was  formed.  Separate  dc  power  supplier  were 
used  to  supply  heater  power  to  the  filament  and  the  bias  voltage  across  the 
plates  number  1  and  number  2.  The  filament  was  maintained  at  ground  otential 
was  plate  number  1. 

Metastables  created  in  the  discharge  source  exit  through  the  0.0.8  cm 
diameter  orifice  that  was  drilled  into  the  chamber  bulkhead  and  enter  a  detec¬ 
tor  chamber  where  they  are  monitored  by  an  Auger  detector.  A  ca'ibrateo 
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Figure  3  Photograph  of  Pumping  System  and  Vacuum  Chamber 
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Figure  4  Schematic  of  Pumping  System  and  Vacuum  Chamber 


Figure  5  Photograph  of  Discharge  Source 
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Figure  6  Schematic  Layout  of  Discharge  Source  and  Detector 


22 


Schulz-Phelps  type  ionization  gauge  was  used  to  monitor  the  total  beam  density 
in  the  detector  chamber.  Oeflector  plates  attached  to  the  chamber  bulkhead 
were  used  to  remove  charged  particles  from  the  beam,  and  a  capillary  tubing 
used  for  injecting  scattering  gas  was  interposed  between  the  deflector  plates 
and  Auger  detector.  Provision  was  also  made  for  a  LiF  beam  flag  which  could 
be  positioned  in  front  of  the  Auger  detector.  Ihe  Auger  detector,  ionization 
gauge  and  gas  cell  were  all  mounted  on  a  translation  stage  that  provided  mo¬ 
tion  perpendicular  to  the  beam  axis  so  that  either  the  Auger  or  the  ionization 
detector  could  be  positioned  on  the  beam  line.  A  picture  of  the  detectors  is 
shown  in  Figure  7. 

A  schematic  diagram  of  the  Auger  detector  and  the  associated  biasing 
circuit  is  shown  in  Figure  8.  The  Auger  plate  is  made  of  pure  molybdenum  with 
a  Q.05  cm  gold  overcoat  to  provide  an  inert  surface.  Metastables  incident  on 
the  gold  surface  will  eject  electrons  if  the  metastable  energy  is  greater  than 
the  work  function  of  gold  (4.3  eV).  The  highly  transparent  attractor  grid  is 
biased  positively  with  respect  to  the  Auger  plate  so  that  the  ejected  Auger 
electrons  are  completely  removed  from  the  plate  resulting  In  a  positive  cur¬ 
rent  monitored  by  a  high  sensitivity  electrometer.  The  guard  plates  insure 
that  leakage  currents  between  the  various  elements  of  the  Auger  stack  are  min¬ 
imized.  This  is  essential  for  the  me„,s«-ement  of  the  low  metastable  currents 
that  were  obtained  in  these  experiments. 

The  Schulz-Phelps  ionization  gauge  calibration  was  verified  by  filling 
the  detector  chamber  with  argon  gas  and  comparing  the  ionization  gauge  mea¬ 
surements  with  that  of  an  absolute  reading  McLeod  gauge.  The  results  of  this 
calibration  are  shown  in  Figure  9.  The  ionization  gauge  measurements  were  ob¬ 
tained  by  scaling  the  ionization  gauge  output  reading  by  the  factor  1/1.29, 
which  accounts  for  the  difference  in  the  ionization  efficiency  of  argon  and 
nitrogen.  The  calibration  points  are  indicated  by  the  solid  dots,  and  the 

solid  line  with  unity  slope  indicates  what  perfect  agreement  would  look  alike. 

-4 

The  ion  gauge  readings  are  reasonably  accurate  down  to  about  10  torr,  the 
lowest  pressure  that  the  McLeod  gauge  can  measure.  Based  upon  the  scatter 
about  the  solid  line,  the  ion  gauge  readings  are  probably  correct  to  better 

than  a  factor  of  two  in  this  pressure  range.  However,  under  the  operating 

-5  -6 

conditions  of  the  present  experiments,  pressures  in  the  10  -  10  torr 
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Figure  7  Photograph  of  Detectors 
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Figure  8  Schematic  of  Auger  Detector  and  Operating  Circuit 
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range  were  measured.  We  assume  that  the  error  limit  obtained  at  the  higher 
pressures  can  be  used  in  this  lower  pressure  regime. 

0.  DIAGNOSTICS  AND  RESULTS 
1.  Discharge  Diagnostics 

The  initial  experiments  involved  measuring  the  discharge  current  en¬ 
hancement  due  to  partial  ionization  of  an  argon  atom  beam.  Figure  10  indi¬ 
cates  the  total  current  drawn  between  the  filament  and  plate  number  2  in  the 
absence  of  any  ionization.  The  data  was  obtained  for  a  0.0127  cm  diameter 
filament  at  an  operating  temperature  of  2138*K.  The  temperature  measurement 
was  obtained  by  optical  pyrometry  including  a  correction  for  transmission 
losses  through  a  quartz  viewing  port.  The  bias  plate  separation  was  1  cm. 

The  space  charge  limited  behavior  of  the  I-V  characteristic  is  clearly  indi¬ 
cated  by  the  slope  of  curve  drawn  through  the  data  points  (crosses).  At  a 
fixed  bias  voltage  (+35  V),  the  total  current  was  then  measured  as  a  function 
of  the  pressure  inside  the  source  gas  nozzle.  This  data  is  displayed  in  Fig¬ 
ure  11.  The  current  rises  only  slowly  for  pressures  up  to  about  150  tor*'  at 
which  point  the  current  increases  dramatically  (by  almost  two  orders  of  magni¬ 
tude)  over  a  narrow  pressure  range,  and  then  rises  slowly  for  higher  source 
pressures.  Visual  inspection  indicates  that  a  discharge  is  ignited  in  this 
high  current  region  as  is  evidenced  by  a  pale  blue  glow  in  the  vicinity  of  the 
gas  nozzle  source.  Similar  results  were  found  for  both  pure  krypton  and  pure 
nitrogen  beams  except  that  the  transition  to  the  high  current  mode  occurred 
for  nozzle  pressures  aruuud  iuu  torr  for  Kr  and  180  torr  for  nitrogen.  No 
discharge  in  helium  was  observed  for  source  pressure  up  to  1000  torr. 

During  measurements  of  the  metastable  production  in  argon,  krypton  and 
nitrogen,  source  pressures  in  excess  of  200  torr  were  not  investigated  due  to 
an  excessive  gas  load  which  could  not  be  handled  by  the  source  chamber  pumpiny 
system.  At  this  pressure  the  gas  flow  from  the  0.018  cm  diameter  nozzle  is  of 
the  order  of  a  few  torr-l/sec. 


2.  Detector  Diagnostics 

The  deflector  plate  voltage  and  the  attractor  grid  voltage  were  opti¬ 
mized  for  metastable  collection  by  monitoring  the  positive  current  emanating 
from  tne  Auger  collector  plate  as  a  function  of  each  of  these  voltages.  With 
the  discharge  lit  and  a  nozzle  pressure  of  about  160  torr  argon,  the  Auger 


CURRENT  (AMPS) 


Figure  10  I-V  Characteristics  of  Discharge  Source  in  the  Absence  of  Ionization 


DISCHARGE  CURRENT  (  AMPS ) 


Figure  11  Discharge  Current  as  a  Function  of  Nozzle  Source  Pressure 
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signal  is  observed  to  be  insensitive  to  deflector  plate  voltages  gre-ate  than 
about  50  V.  A  similar  result  can  be  obtained  by  biasing  the  whole  Auger  stack 
5U  V  negative  with  respect  to  the  grounded  discharge  source  suggesting  that 
high  energy  electrons  are  the  primary  charged  species  carried  in  the  beam.  As 
is  indicated  in  Figure  12,  an  attractor  grid  voltage  of  about  +  20  V  with  re¬ 
spect  to  the  Auger  plate  is  required  to  insure  that  all  the  ejected  Auger 
electro..  ',re  removed  from  the  Auger  plate.  Higher  attractor  grid  voltages 
are  avoided  to  insure  that  the  hign  energy  discharge  electrons  do  not  impinge 
on  the  Auger  plate  and  create  secondary  electrons. 

3.  Results 

Using  these  optimal  bias  voltages,  the  Auger  signal  was  measured  as  a 
function  of  the  source  nozzle  pressure.  Ihese  results  are  displayed  in  Figure 
13  for  conditions  identical  to  Figure  11.  A  comparison  of  these  two  figures 
indicates  that  the  Auger  current  behavior  closely  follows  that  of  the  source 
discharge  current.  Several  additional  checks  were  made  to  insure  that  the 
Auger  current  measured  under  “discharge  ori"  conditions  is  a  true  measure  of 
the  metastable  beam  flux  and  is  not  contaminated  by  contributions  from  dis¬ 
charge  photons  or  charged  particles. 

The  photon  contribution  from  a  Kr  discharge  was  estimated  be  observing 
the  attentuation  of  the  Auger  signal  when  a  LiF  beam  flag  is  interposed  in  the 
beam  line  to  block  all  particles  from  hitting  the  detector.  LiF  crystal  effi¬ 
ciently  transmits  over  a  very  wide  optical  wavelength  range  (10  percent  cutoff 
points:  100  and  9000  nm)  including  the  krypton  resonance  line  (*Pj  >  *SQ) 
at  1165  A  where  the  LiF  optical  transmission  is  about  20  percent.  Using  a 
krypton  discharge,  the  Auger  signal  attenuates  by  more  than  a  factor  of  100 
when  the  beam  flag  is  inserted  directly  in  front  of  the  Auger  detector.  This 
suggests  that  the  photoemission  contribution  to  the  Auger  signal  is  at  most  a 
few  percent  for  the  krypton  metastable  measurements.  The  similarity  in  the 
operating  characteristics  of  the  krypton,  argon  and  nitrogen  discharges  sug- 
gets  that  photoemission  also  contributes  negligibly  to  the  argon  and  nitrogen 
Auger  signals. 

The  presence  of  argon  metastables  in  the  beam  was  uniquely  identified  by 
a  Penning  ionization  technique.  The  lowest  lying  argon  metastable  (3 P ^  at 
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AUGER  CURRENT  (AMPS) 


Figure  13  Auger  Current  as  a  Function  of  Nozzle  Source  Pressure 
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11.5  eV)  is  sufficiently  energetic  to  Penning  ionize  C2H2  (ionization  po¬ 
tential  «  11.4  eV)  but  is  insufficient  to  Penning  ionize  09  (ionization  po- 
tential  ■  12.06  eV).  On  the  other  hand,  Ar  (15.76  eV)  can  rapidly  charge 
transfer  to  0.,.  A  comparison  of  the  ion  currents  collected  when  either 
CgHg  or  02  are  added  to  the  metastable  beam  should  give  a  relative  indi¬ 
cation  of  the  Ar  and  Ar  (metastable)  beam  content.  This  measurement  was 
accomplished  by  introducing  a  tiny  flow  of  either  gas  through  a  capillary  tube 
into  the  region  between  the  deflector  plates  located  directly  ahead  of  the 
Auger  plate  (see  Figure  7).  The  flow  was  directed  nominally  perpendicular  to 
the  beam  line,  and  the  resultant  ion  currents  were  collected  on  the  deflector 
plates  and  measured  with  a  sensitive  electrometer.  As  CgH^  is  added,  the 
current  to  the  deflector  plates  increases  by  about  a  factor  of  30  at  which 
point  the  current  begins  to  saturate.  Under  similar  conditions,  the  Auger 
current  is  reduced  by  a  factor  of  about  ten  due  to  scattering  and  Penning 
ionization.  For  comparison,  a  negligible  increase  in  the  deflector  plate 
current  is  observed  when  a  comparable  amount  of  oxygen  is  added  sufficient  to 
reduce  the  current  to  the  Auger  plate  by  a  similar  factor  of  ten.  These  re¬ 
sults  indicate  that  the  ion  content  of  the  beam  is  negligible  compared  to 
metastable  content,  and  that  some  fraction  of  the  metastables  is  comprised  of 
Ar  (  P^).  We  were  unable  to  perform  similar  tests  on  krypton  and  nitrogen 
metastable  beams.  A  potential  problem  for  the  N2  metastable  measurements  is 
our  current  inability  to  distinguish  between  atomic  and  molecular  metastable 
states  in  nitrogen.  The  atomic  metastable  states  could  be  formed  by  dissocia¬ 
tive  excitation  processes  in  N2  discharges.  Although  the  low  lying  N(4P) 
and  N(4D)  states  have  insufficient  energy  to  e.iect  Auger  electrons  from  the 
gold  surface  (recommended  work  function  =  4.3  eV),  higher  energy  N  atom  meta¬ 
stable  levels  could  contribute.  The  relative  contribution  of  these  states  to 
our  measured  Auger  current  is  not  known  at  present. 

Having  concluded  that  the  argon  Auger  current  is  a  true  measure  of  the 
argon  metastable  beam  flux,  the  discharge  source  parameters  were  varied  in 
order  to  optimize  the  metastable  flux.  The  discharge  bias  voltage,  discharge 
current  and  the  nozzle  source  pressure  are  the  parameters  of  primary  interest 
here.  Figure  11  has  already  indicated  that  there  is  only  a  slight  increase  in 
metastable  flux  for  nozzle  pressures  in  excess  of  that  required  to  ignite  the 
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discharge.  Therefore,  in  light  of  this  observation  as  well  as  for  the  need  to 
remain  within  the  pumping  speed  limitations  of  the  vacuum  system,  we  have  held 
the  nozzle  pressure  constant  at  -  ISO  torr  (for  argon),  -  ICO  torr  (for  kryp¬ 
ton),  -  180  torr  (for  nitrogen)  in  these  studies. 

At  fixed  filament  power,  the  discharge  voltages  at  which  the  metastable 
signal  optimizes  are  found  to  be:  30  to  40  V  for  krypton  or  argon  and  50  to 
60  V  for  nitrogen.  Representative  data  for  the  optimization  in  a  nitrogen 
beam  are  given  in  Figure  14.  Ihis  variation  is  similar  to  that  obtained  in 
argon  and  krypton  and  is  due  to  both  variations  in  both  the  electron  energy 

and  electron  current.  Ihese  effects  can  be  decoupled  by  varying  the  filament 

power  at  fixed  discharge  voltage.  This  data  is  shown  in  Figure  IS  for  both 
argon  and  krypton  and  nitrogen  beams.  In  contrast  to  the  data  in  Figure  14, 
the  metastable  signal  varies  monatonically  with  discharge  current  for  each  gas 
studied.  Ihe  Auger  current  appears  to  be  limited  only  by  the  filament  burnout 
temperature  (-  2500*K)  which  is  attained  at  about  30  W  (Vp  -  7  V,  Ip  «■ 

4.3  amps)  heater  power  corresponding  to  the  high  end  of  the  abscissa  in  Figure 

15.  A  comparison  of  Figure  14  and  Figure  15  then  suggests  that  the  rollover 
at  high  discharge  voltage  can  be  directly  attributed  to  an  electron  energy  ef¬ 
fect.  We  speculate  that  for  discharge  voltages  below  the  peak  in  Figure  14, 
the  increase  in  the  ionization  cross  section  with  increasing  electron  energy 
dominates  over  any  fall-off  in  the  metastable  excitation  cross  section  leading 
to  a  net  increase  in  the  metastable  signal.  As  the  discharge  voltage  in¬ 
creases  beyond  the  peak,  the  decrease  in  the  metascable  cross  section  begins 
to  dominate  over  the  slower  rising  ionization  cross  section  leading  to  a  net 
decrease  in  the  metastable  signal. 

Several  quantities  are  required  to  convert  the  Auger  current  measure¬ 
ments  in  Figure  15  into  a  metastable  density  value  at  the  exit  of  the  dis¬ 
charge  zone:  the  secondary  emission  coefficient,  n,  for  metastable  impact  on 
a  gold-coated  surface,  the  velocity  distribution  of  the  metastable  beam,  the 
angular  distribution  of  metastables  exiting  from  the  discharge  and  the  solid 
angle  subtended  by  the  detector.  Of  these  only  the  Ar  (  P_  )  secondary 

c90 

emission  coefficient  and  the  detector  solid  angle  are  known  for  these  experi¬ 
ments.  The  secondary  emission  coefficient  has  been  reported  to  be  0.66  ±0.1 
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Figure  14  Auger  Current  as  a  Function  of  Discharge  Voltage 
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Figure  15  Kr,  Ar  and  Ng  Metastable  Signals  as  a  Function  of  Discharge 
Current 
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for  Ar  ( 3Pg  o)^50^  The  detector  solid  angle,  &  ,  was  established  by 
locating  a  small  aperture  (0.079  cm3)  directly  'n  front  of  the  Auger  plate 
at  a  distance  of  8  cm  from  the  discharge  exit  orifice.  The  angular  and 
velocity  distributions  of  the  metastable  beam  need  to  be  obtained  experi¬ 
mentally  due  to  the  possibility  of  hydrodynamic  flow  at  the  discharge  exit 
orifice.  However,  these  measurements  are  difficult  and  were  not  attempted 
here. 

We  take  what  information  is  available  to  us  and  report  the  metastable 
flux  per  unit  solid  angle  (a  quantity  analogous  to  brightness  in  optics),  F, 
given  by  the  expression 


!W 

&»Asne 


(6) 


M 

where  I  is  the  measured  Auger  current  with  the  aperture  in  place,  A$  is 
the  discharge  exit  orifice  cross-sectional  area  and  e  is  the  electron  charge. 
We  find  that  the  IM  for  argon  is  a  linear  function  of  the  discharge  source 
current,  Irf,  given  by 


IM  «  1.4  x  10“9  Id  (in  amps) 


(7) 


for  Id  <  0.5  A.  Using  this  value,  the  quantity,  F,  was  evaluated  for  as;  = 

1.2  x  10“3  sr,  A$  ■  1.14  x  10“3  cm3  and  rj  -  0.66  to  be  about  1016  meta- 
stables/cm^-  sr  for  an  argon  discharge.  Using  this  value  and  our  measure¬ 
ments  of  the  atomic  beam  density,  we  estimate  a  fractional  beam  content 
of  the  order  of  0.1  percent.  Once  the  secondary  emission  coefficients  for  N? 
and  Kr  are  known,  the  metastable  yields  for  these  gases  can  be  scaled  from  F 
in  argon  using  the  data  in  Figure  15. 

E.  RECAP  OF  SOURCE  AND  DE1EC10R  OPERATING  CONDITIONS 

Ihe  optimal  operating  conditions  for  the  discharge  source  and  Auger  de¬ 
tector  are  reviewed  below: 


Uin  Dunning,  F.B.,  and  Smith,  A.C.H.,  J.  Phys.  84,  1969  (1971). 
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a)  Detector  Conditions 

Attractor  Bias  »  +20  V  relative  to  Auger  plate 
Deflector  Bias  >  50  V 

Auger  Plate  =  -50  V  relative  to  case  ground 
Auger  stack  bias  ■  -50  V  relative  to  case  ground 

b)  Source  Conditions 

Nozzle  Pressure  =  100  torr  (Kr) ;  150  torr  (Ar);  180  torr  (Ng) 
Discharge  Voltage  a  +30-40  V  (Kr  and  Ar);  +50-60  V  (N^) 

Discharge  Current  <  700  mA 

Filament  lemperature  (0.0127  cm  diameter)  <  2500*K 
F.  PROPOSAL  FOR  MEASUREMENT  OF  ABSOLUTE  METASTABLE  DENSITY 

We  propose  using  the  Penning  ionization  technique  described  previously 
to  make  absolute  density  measurements  of  argon  metastables.  The  technique  re¬ 
quires  making  absolute  measurements  of  the  electron  currents  generated  from 
the  process: 


Ar*  +  M  »  M+  +  e  +  Ar 


(8) 


where  Ar*  is  a  metastable  excited  state  and  M  is  an  appropriate  collision 
partner  that  can  be  Penning  ionized.  For  detecting  Ar  metastables,  the  ioni¬ 
zation  potential  of  this  species  must  be  £  11.55  eV,  the  energy  of  the  lowest 
Ar  metastable  state.  The  electron  current  generated  via  reaction  (8)  will 
tend  to  a  saturated  value,  1^,  as  the  M  gas  density  is  raised.  At  this 
point  all  the  Ar*  metastable  species  will  be  ionized.  The  value  of  Ifi 
can  then  be  related  to  the  number  density, pm,  via: 


PM 


(9) 


where  V£  is  the  Penning  electron  velocity,  A  is  the  electron  collector  area 
and  e  is  the  electronic  charge.  In  effect,  this  procedure  converts  one  spe¬ 
cies  (metastables)  whose  density  is  difficult  to  directly  measu.'e  into  another 
species  (electrons)  the  density  of  which  is  relatively  easy  to  measure. 
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The  metastable  density  can  be  obtained  from  measurements  of  the  total 
current  and  the  electron  energy  distribution.  This  latter  quantity  can  be  ob¬ 
tained  using  an  electrostatic  energy  analyzer.  An  accurate  measurement  of  1* 
would  require  a  Penning  electron  collection  efficiency  approaching  100  nercent. 
Electron  scattering  into  angles  outside  the  field  of  view  of  the  electron  col¬ 
lection  optics  would  be  the  major  effect  that  could  reduce  the  efficiency  be¬ 
low  unity.  The  collection  efficiency  would  be  maximized  by  using  an  M-gas  for 
which  the  cross-section  ratio  Op/o^,  the  ratio  of  the  Penning  cross  section 
to  the  momentum  transfer  cross  section,  is  large.  This  would  insure  minimal 
amount  of  e-M  lar/e  angle  scattering.  The  collection  efficiency  could  also  be 
optimized  by  choosing  the  M-gas  such  that  the  product  electron  would  have 
little  excess  translational  energy. 

C9Ho  is  a  reasonable  M  gas  candidate.  The  energy  difference  between 
3  ^  £  + 

Ar(  P2)  and  C2H2  is  only  0.15  eV.  The  momentum  transfer  (and  rotational 

excitation)  cross  sections  are  not  known  although  it  is  expected  that  oM(C2H2) 

is  small  due  to  its  zero  dipole  moment  and  small  polarizability.  On  the  other 

hand,  the  Ar(  Po)/C,H0  Penning  ionization  cross  section  has  been  measured 
-15  6 

to  be  >  10  cm  at  high  energy  (>  100  eV),  and  the  cross  section  is  ris- 

(51)  “ 

ing  towards  lower  energy.'  ' 


51.  Moseley,  O.T.,  Peterson,  J.R.*  Lorents,  D.C.,  and  Hollstein,  M.,  Phys. 
Rev.  A6,  1025  (1972). 
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IV.  HE THODS  FOR  PRODUCING  HIGH  THROUGHPUT  SOURCES  OF  HE TAS TABLE 

ATOMS  AND  MOLECULES 


The  objective  of  this  phase  of  the  program  is  to  investigate  the  possi¬ 
bility  of  developing  sources  that  can  generate  high  flux  levels  of  metastablo 

atoms  and/or  molecules.  An  arbitrary  preliminary  design  goal  of  the  order  of 
PI  22  -2-1 

lu  -  10  metastables  cm  sec  is  assumed.  An  additional  constraint 
requiring  flow  velocities  that  do  not  exceed  the  speed  of  sound  insures  that 
the  internal  energy  of  the  metastables  greatly  exceeds  the  translational  ki- 
netic  energy.  At  a  typical  sound  speed  of  the  order  of  10  cm/ sec,  the  re¬ 
quired  metastable  density  levels  are  of  the  order  of  lO1^  -  1018  em“3,  a 
range  that  is  orders  of  magnitude  larger  than  most  any  metastable  source  (ex¬ 
cept  one)  with  which  we  are  familiar.  The  one  exception  is  the  ) 
chemical  metastable  generator  that  is  used  in  the  0o(1a„)/I  chemical 

1  *■  i  ■) 

laser  system.  ag)  densities  of  the  order  of  **  10A/  cm-0  have 
been  reported. 

Four  different  approaches  were  considered  in  this  study.  Discharge  pro¬ 
duction  of  the  metastable  Ng( A)  state  was  investigated  using  an  e-beam  sus¬ 
tained  discharge  to  generate  volumetric  pumping  at  the  high  energy  loadings 
needed  to  provide  a  high  flux  level.  The  ability  to  independently  vary  the 
e-beam  rurrent  and  discharge  voltage  to  optimize  the  A-state  production  is  a 
decided  advantage  of  the  e-beam-sustained  discharge  over  more  conventional 
self-sustained  discharge  sources.  On  the  other  hand,  the  issue  of  the  hot 
electrons,  which  both  produce  the  A-state  anr,  unfortunately  destroy  it  by 
cumulative  ionization  and  dissociation  processes,  arises  in  all  discharge 
pumping  techniques  and  cannot  be  avoided.  An  alternative  approach  which 
attempts  to  circumvent  the  hot  electron  problem  is  proauction  of  metastables 
via  electron-ion  recombination  processes  which  are  most  efficient  for  cold 
electrons  due  to  the  inverse  electron  temperature  dependence  of  most  electron- 
ion  recombination  rate  constants.  We  have  investigated  this  concept  using  a 
magnetic  annular  arc  (MAARC)  generator  as  a  high-temperature  thermal  ioniza¬ 
tion  source  of  ions  and  electrons.  The  recombination  of  ions  and  electrons  is 
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allowed  to  take  place  in  a  cool  region  that  is  spatially  removed  from  the  ion 
source.  The  plasma  is  physically  transported  from  hot  to  cold  regions  by  a 
flow  and  confining  magnetic  fielu 

Experience  in  atomic  and  molecular  physics  suggests  that  optical  excita¬ 
tion  techniques  should  be  expected  to  be  more  efficient  and  selective  than 
processes  using  electrons.  Generally  speaking,  the  optical  selection  rules 
and  photon  energy  resonance  requirements  lead  to  preferential  excitation  of 
fewer  levels  and  decidely  fewer  loss  processes  than  occurs  for  nonresonant 
electron  impact  excitation.  As  a  result,  we  have  investigated  two  purely  op¬ 
tical  metastable  production  techniques.  The  first  uses  a  two-color  optical 
pumping  scheme  for  production  of  low  lying  metastable  levels  in  the  heavier 
rare  gases  (krypton  and  xenon).  The  second  uses  a  stimulated  Raman  technique 
to  produce  vibrational ly  excited  Hg  and  Ng. 
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A.  STIKULATtiO  RAMAN  SCATTERING  PRODUCTION  OF  «2  (v«l)  find  N?  (v«.l) 

1.  Introduction  and  Background 

Stimulated  Raman  scattering  (SkS)  involves  an  inelastic  scattering 
mechanism  in  which  a  photon  incident  on  an  atom  or  molecule  causes  that  atom 
or  molecule  to  undergo  a  change  in  internal  energy  and  emit  a  photon  that  is 
shifted  in  wavelength  from  the  incident  photon  by  that  exact  change  in 
internal  energy.  In  the  case  of  an  atom,  the  internal  energy  is  in  the  form 
of  an  electronic  transition  (see  Figure  16).  For  a  molecule,  the  change  in 
internal  energy  can  be  vibrational  or  rotational.  We  intend  to  make  use  of 
SRS  to  produce  significant  fractional  populations  of  vibrational ly  excited 
diatomic  molecules  at  high  densities  to  meet  the  de  ired  hiqh  flux 
requirements. 

The  molecular  Raman  process,  involving  a  very  non-resonanf  scatterinr, 
process,  has  a  smaller  Raman  cross  section  than  is  typical  for  the  atomic  case 
anti,  therefore,  experiments  are  typically  carried  out  at  high  pressure.  In 
fact,  pressure  is  one  of  the  most  significant  parameters  affecting  the 
molecular  SRS  process.  This  is  because  the  variation  in  pressure  will  affect 
both  the  number  density  of  the  scatterer  and  the  spontaneous  Ramon  linewidth. 
Since  the  gain  is  equal  to  the  Raman  cross  section  times  the  density  tines  the 
length,  i.e., 

Gain  n  chi 

and  the  Raman  cross  section  is  inversely  proportional  to  the  Raman  linewidth 

for  the  case  of  forward  SRS,  a  consideration  of  the  pressure  influence  on  the 

gain  is  important.  This  will  enable  us  to  evaluate  the  quantitative 

performance  of  this  Raman  approach  for  a  variety  of  system  total  pressures. 

For  Hj,,  there  are  three  distinct  pressure  regimes  for  the  Raman 

linewidth  (see  Figure  17).  At  very  low  pressures,  the  linewidth  is  primarily 

—  2 

determined  by  Doppler  broadening  avq  c  2v/C  >'2  1  n2  kT/M  =  3.6  x  10 
cm"1  for  forward  scattering.  In  the  collision  narrowino  regime, 

=  4*D/x  1/C  =  5.74  x  10~3  cm_1/N(Amg)  where  D  is  the  diffusion 
coefficient  of  Hg  in  (D  =  1.35),  C  is  the  speed  of  light  and 
avcn=  avq  at  N  =  0.27  Amg  or  205  torr.^52^  At  slightly  higher  pressures,  the 

52.  Hurray,  J.R.  and  Javan,  A.,  0.  Molec.  Spectry  42,  1  (1972);  Dicke,  R.H., 
Phys.  Rev.  89,  472  (1953).  ~ 
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Figure  16  Stimulated  Raman  Scattering  (SRS) 
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Figure  17  Pressure  Effect  of  H2  Raman  Linewidth 
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forward  scattering  Raman  linewidth  decreases  due  to  collision  narrowing 
effects  resulting  from  the  frequent  velocity-changing,  but  not  phase-changing, 
collisions.  In  this  regime,  the  linewidth  becomes  inversely  proportional  to 
the  density  of  the  collision  partner  (i.e.,  Hg)  and,  therefore,  the  gain 
scales  as  the  density  squared.  At  still  higher  pressures  where  A\>CN  =  AvpB 
i.e.,  the  linewidth  becomes  pressure  broadened,  the  gain  becomes  independent 
of  pressure. 

We  can  how  calculate  the  gain  in  each  of  these  three  regimes,  first 
defining  a  g(x)  such  that 


g(x) 


f(x) 


Av 


with  a w  independent  of  pump  wavelength.  Therefore, 

9PBw  - 

AVod 


Av 


PB 

9cn(*)  ■  9PB(*)  -r—  =  9pg(x) 
AvCN 


1.5  x  10~3  N 

9774  x  10“3/N 


»  0.154  N2  gpg(x) 


g0(x)  «  gPB(x)  9PB(*)  [3*5  x  W  =  °'0417N  W*) 


The  value  of  gDO(x)  in  cm/Vf  as  a  function  of  wavelenqth  is  plotted  in 

(531 

Figure  18  normalized  to  the  data  of  Bloembergen  et  al.  '  The  pressure 
braodened  linewidth  of  1.5  x  10"3N  was  taken  from  Hurray  et  al.^4^  Using 
these  relationships,  with  various  laser  pump  sources,  we  arrive  at  the 
following  values: 


SD 

gCN 

gPB 


Ruby  (x  =694  nm) 

6.3  x  HT11  N 

-10  2 

2.3  x  1C  iU  r 


KrF  (x  =  248  r.m) 
3.6  a  10"10  N 
1.3  x  10"9  U2 


ArF  (x  =  193  nm) 
6.7  x  10"16  N 
2.5  x  10“9  N2 


1.5  x  10 


-9 


8.6  x  10 


-9 


1.6  x  10 


-8 


53.  Bloembergen,  N.,  Bret,  G.,  Lallemand,  P.,  Pine,  A.  and  Simova,  P.,  IEEE 
0.  Quant.  Electron  QE^,  197  (1967). 

54.  Murray,  J.R.,  Goldhar.  J.,  Eimer,  D.,  and  Szoke,  A.,  IEEE  0.  Quant. 
Electron  QE-15,  342  (1979). 
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Figure  18  First  Stokes  Gain  for  Forward  SRS  in  H2 
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Similarly,  one  can  explore  the  analogous  processes  for  SRS  in  nitrogen 
but  unfortunattiy  little  information  is  available  with  regard  to  the  effect  of 
pressure  on  the  Raman  linewidth  in  the  regime  of  interest.  Murray  et 
al.^^  suggest  that  for  pressures  of  a  few  atmospheres  or  more,  the 
linewidth  is  probably  independent  of  density.  At  very  high  pressures  (>  200 
Amg)  the  linewidth  begins  to  narrow  such  that  at  liquid  densities  it  is  only 
about  0,06?  erf1.  For  this  reason,  most  SRS  studies  involving  nitrogen  have 
been  carried  out  using  liquid  nitrogen,  which  is  unacceptable  for  the  present 
application.  If  the  Raman  linewidth  is  indeed  independent  of  density  for 
pressures  below  1  atm,  then  the  gain  scales  with  density.  Using  the  published 
value  for  g  at  488  nm  and  scaling  with  wavelength  to  x  »  560  nm  (see  below  for 
choice  of  appropriate  wavelength)  we  obtain 

g  *  i .38  x  10"12  N, 

which  is  considerably  smaller  than  those  for  Hg  (see  above). 

These  cross  sections  can  now  be  used  to  calculate  the  expected  number 
densities  of  vibrationally  excited  Hg  or  that  can  be  produced  by 
co-propagating  (i.e.,  coherently  mixing)  both  pump  and  Stokes  photons,  which 
act  to  shift  population  from  the  lowest  v  a  0  state  (i.e.,  N)  to  the  v  n  1 
state  (i.e.,  N*) .  At  this  time,  let  us  assume  the  Stokes  photons  are  produced 
arbitrarily.  A  rate  expression  describing  this  excitation  process  is 

dN*  X  M 

d r=  °srs  fis N 


where  oSRS  is  the  SRS  cross  section  in  enr  and  is  the  flux  of  photons 
at  the  Stokes  frequency.  Assuming  no  saturation  effects  (i.e.,  small  changes 


1n  N'  ‘pump  and  ‘s'- 


>S  N  T pulse 


Substituting 

°SRS  E  °'SRS  ]P» 

where  o'^  is  the  SRS  "cross  section"  in  crn^/W,  we  obtain 

N*  =  °SRS'  1 P  ^S  N  T pulse 
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Now  c^'  N  n  g(cm/W),  which  we  have  tabulated  above,  and  $s  n  ^ s^vs» 
therefore 


and  values  of 

g0/hv$ 

9CN/hvs 

9pi^hus 


g/hv.  can  be  tabulated  as  before.  For  H«, 


Ruby  (x  a  694  nm) 

3.2  x  108  N 

1.2  x  109  N2 
7.8  x  109 


KrF(x  ■  248  r.m) 
5.0  x  108  N 
1.8  x  109  N2 
1.2  x  1010 


ArF  (x  a  193  nm) 

7.0  x  108  H 
q  2 

2.6  x  10'  N* 

1.7  x  1010 


5 

whereas  g/hv  for  nitrogen  is  4.5  x  10  N  at  x„.  »  560  nm.  From  these 

3  s  pump 

data,  it  is  apparent  that  the  quantity  g/hvs  is  very  weakly  wavelength 

dependent.  This  arises  because  g  ~  vsg'(x),  and  g'(x)  has  a  wavelength 

dependence  which  is  only  due  to  the  resonance  denominator,  which  is  very 

non-resonant  for  all  of  these  pu:..ps  and,  therefore,  a  very  weak  function. 

Although  ArF  has  minor  gain  advantage,  its  short  wavelength  poses  optics 

damage  considerations  and,  therefore,  we  shall  pursue  our  calculations  using 

KrF  (x  «*  248  nm)  as  the  pump  for  H2.  We  will  also  assume  we  will  generate 

the  needed  Stokes  radiation  (x  =  277  nm)  in  a  high-pressure  hydrogen  cell. 

2.  Proposed  Experimental  Set-Up  and  Predicted  Yields 

Based  on  experience  at  AERL  with  commercially  available  KrF  lasers  and 
(55) 

SRS  conversion  in  Hg,  the  following  values  for  the  needed  parameters 
are  considered  to  be  realistic  estimates.  For  an  unstable  cavity  operated 
Lumonics  (Model  TE-261-2)  laser,  energy  expected  at  KrF  wavelengths  is  80  mO 
in  a  8  nsec  pulse.  From  our  SRS  work  8 ^  and  that  of  others/88^  we 
project  we  can  easily  convert  50  percent  of  this  energy  into  a  number  of 
Stokes  and  anti-Stokes  lines  with  most  (75  percent)  going  into  the  desired 
S^  radiation.  Therefore,  after  the  pump  has  been  converted  in  t!*e  hioh 


55.  Trainor,  D.W.,  Hyman,  H.A.,  Itzkan,  I.  and  Heinrichs,  R.M.,  App.  Phys. 
Lett,  37,  440  (1980). 

56.  Loree,  T.R.,  Sze,  R.C.  and  Barker,  D.L.,  App.  Phys.  Lett.  31,  37  (1977); 
Loree,  T.R.,  Sze,  R.C.,  Barker,  D.L.  and  Scott,  P.B.,  IEEE^J.  Quant. 
Electron  QE-15,  337  (1979). 
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pressure  cell  (see  Figure  19),  we  have  Ep  *  <10  f»\)  and  Es  ■  30  mJ.  This 
30  mO  of  277  nm  radiation  corresponds  to  <1.2  x  1016  photons,  thereby 
establishing  the  maximum  number  of  excited  Hg(v«l)  states  that  can  be 
achieved  with  this  pump.  This  number  and  the  desire  to  have  high  fractional 
excitation  places  constraints  on  operating  conditions. 

Since  Ip  »  Ep/TpA  and  I$  ■  E$/*pA,  where  A  is  the 
cross-sectional  area  of  the  beam,  then 

"*  ■  to  &)  (*?)  • 


Also,  with  A  ■  isw“,  then  N*  for  Hg  for  each  pressure  regime  can  be 
tabulated  as  a  function  of  beam  waist  in  the  farfield  (i.e.,  focal  spot 
size).  For  these  assumed  experimental  conditions,  we  have 


2.75  x  10“7  II 


N  <  7.3  x  1018  cm"3 


..  .  3.75  x  iCf26  .,2 

V8 — ^ —  N 


7.2  x  1018  <  N  <  7  x  1019  cm"3 


1.8  x  10' 


N  >  7  x  1019  cm"3 


where  the  beam  waist  (for  a  diffraction  limited  beam)  is  given  by 

..  1.22  f  X 


where  f  is  the  fecal  length  of  the  optics  employed,  x,  is  the  wavelength  and  0 
is  the  diamet  of  the  incident  laser  beam.  The  depth-of -field  over  which  the 
beam  is  within  a  factor  of  the  Si  of  this  beam  waist  is  given  by 

„w2 

zv  =  -j—  (i.e.,  the  Rayleigh  Range) 

Therefore,  the  irradiated  volume  is  approximately 

\i  i  2  nV* 

Vol  =  JlW  z  *  — — 

V  A 
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Figure  19  Experimental  Setup  for  H2  SRS 
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Tne  number  of  excited  states  in  this  volume  (for  the  Duppler  regime)  is 


Vol  *  density 


A4  x  2.75  x  1(T7 
__  -  r< 

*  w 


*  excited  states  «  C.098  *  N 

which  must  be  smaller  than  the  number  of  photons  available,  i.e., 
A  AAA  ■*  (I  E  30  x  10“3J —  b  4.2  x  1016 


0.098  *  N  < 


1.6  x  10”1**  4.48 


Assuming  one  is  restricted  to  converting  ~  1/2  the  number  of  photons 
available,  then  H*  Vol  ~  2  x  10^  or  N  <  2  x  10^7  cm  3,  or 

N*  Vol  -  N*  «  2  x  1016 


N*  w4  «  2  x  10i6 

V 


log  N*  +  4  log  w  »  10.75 

A  plot  of  this  equation  is  shown  in  Figure  20.  It  shows  for  10  percent 
excitation  into  the  v  *»  1  level,  the  beam  waist  is  w  •  0.41  mm  »  0.041  cm. 
This  corresponds  to  a  Rayleigh  range  of  190  cm  and  an  irradiated  volume  of 
-  1  cm  .  These  focal  parameters  can  be  accomplished  through  a  system  with 
an  effective  focal  length  given  by 

r  W  D  ~  __  n  1 


1200  cm  for  0  a  1  cm. 


This  long  focal  length  can  be  achieved  in  a  compact  manner  throuqh  use  of  a 
Galilean  type  telescope  arrangement,  which  avoids  breakdown  concerns,  which 
may  occur  with  use  of  a  Keplerian  type. 

We  have  now  defined  a  typical  set  of  experimental  conditions  for  the  low 
pressure  cell.  They  are  summarized  as  follows: 

H,  pressure  =  6.2  torr  =  2  x  10*7  cm"3  for  H9(v=l)/H5  =  10 

c  c  16''  -3 

percent  excitation  implies  Hg  (v=l)  =  0.62  torr  =  2  x  10  cm  for  the 

case  of  a  beam  diameter  »  2w  =  0.08  cm  =  0.8  mm.  Assuming  a  flow  velocity 

which  is  sonic,  v  =  ^KT/M  =  1.5  x  105  cm/sec.  For  this  velocity  and  this 

21  -2  -1 

density,  one  would  anticipate  a  flux  of  ~  3  x  10  cm  sec  ,  which  is 
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near  the  Task  requirements.  In  practice,  both  the  velocity  and  density  at  the 
exit  plane  will  be  reduced  somewhat  but  the  development  presented  here  would 
not  change  significantly.  (Significant  increases  in  flux  can  be  realized 
through  additional  laser  pump  energy  or  different  focal  parameters  (see 
Figure  20)). 

The  Raman  frequency  shift  due  to  the  high-pressure  gas  (K^  «  10  atm) 
in  the  high-pressure  cell  will  te.id  to  drive  the  frequency  difference  between 
the  pump  and  the  Stokes  photon  off  the  resonance  value  for  coherent  excitation 
of  the  low  pressure  Hj,  (p  -  6  torr).  However,  a  full  compensation  can  be 
accomplished  through  the  addition  of  2.73  atm  of  Helium  to  the  high  pressure 
cell  (see  Ref.  57). 

The  ease  of  accomplishing  this  correction  in  Hg  is  to  be  contrasted 
with  Ng,  where  the  gas  phase  vibrational  Raman  shift  has  been  measured  at 
2327  cm-1,  whereas  at  liquid  densities  (a  convenient  high  pressure  cell)  the 
shift  is  2331  cm-*.  This  corresponds  to  a  change  in  frequency  of  4  cm-*, 
which  is  much  too  large  to  be  conveniently  circumvented.  For  this  reason,  we 
have  explored  the  use  of  two  high  power  dye  lasers  to  generate  the  pump  and 
"Stokes"  photons  for  the  calculation  of  SRS  production  of  N2(v«l). 

Current  state  of  '.he  art  commercial  dye  laser  devices  have  demonstrated 
200  mJ  in  10  nsec  at  575  nm  using  a  doubled  yag  laser  to  pump  rhodamln  6  G. 
Similarly,  one  can  expect  near  60  rod  in  10  nsec  at  644  nm  through  use  of 
Rhodamin  640  dye.  These  commercially  available  lasers  would,  thereby,  produce 
the  following  excitation  in  nitrogen. 

2  x  10"9 


As  before,  the  number  of  excited  states  in  a  focal  volume 

Vol  x  density  =  ~ x,*°--  N 

*  vr 

#  excited  states  =  3.07  x  10-4  N 
which  must  be  smaller  than  the  number  of  photons 


3.07  x  10-4  N  <  = 


1.6  x  10-i?*  1.9 


*  2.0  x  10J 


57.  Lallemand,  P.  and  Simova,  P.,  <3.  Molec.  Spectro.  26,  262  (1968). 
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Again,  assuming  one  can  convert  1/2  of  the  photons  available,  then 
tt*  Vol  -  1.0  x  10'  or  N  <  3.2  x  1020  cm”3, 


Therefore, 


U*  Vol 


2  4  ,, 

U*  ~  -  1  x  1017 


or 


log  N*  +  4  Log  w  -  11.8 

This  equation  is  also  plotted  on  Figure  20.  It  shows  that  to  achieve  10 
pircent  excitation  into  v  ■  1,  the  beam  waist  must  be  siqnif icantly  tighter 
than  for  the  HL  case,  i.e.,  w  »  0.12  mm  -  0.012  cm.  This  corresponds  to  a 

c  _3  *» 

Rayleigh  range  of  ~  7  cm  and  an  irradiated  volume  of  -  3  x  10  cm’ . 

Conversely,  the  orifice  oiameter  appropriate  for  10  percent  excitation 

in  H?  will  achieve  only  -  0.7  percent  excitation  in  N?.  For  this  regime, 

the  effective  focal  length  would  be  near  150  cm  f.l.  The  tight  focusino 

10  9 

condition  would  correspond  to  intensities  of  near  4  x  10  W/cm  which  may 
prove  experimentally  troublesome.  For  the  conditions  used  in  the  Hp 
example,  we  have  N?  pressure  B  13  atm.  For  a  beam  diameter  of  0.8  run, 

=  0.7  percent  or  ^(vsl)  =  69  torr.  Assuming  a  flow  velocity 
which  is  sonic  (v  =  4.2  x  10^  cm/sec)  and  these  experimental  conditions, 
one  would  anticipate  a  flux  of  9.4  x  1022  cm"2  sec"1  which  significantly 
exceeds  the  task  requirements. 

We  will  now  consider  the  various  experimental  system  constraints,  having 
demonstrated,  at  least  in  principle,  that  we  can  generate  the  needed  densities 
and  flux  at  high  fractional  excitation  ratios  using  realistic  laser  pumps  for 
both  Hg  and  Np. 

3.  System  Parameters  (Using  Ho  as  an  Example) 

Since  these  lasers  have  pulse  lengths  of  near  10  nsec  the  resultinq  flux 
of  vibrational ly  excited  Hp  exiting  from  the  0.8  mm  orifice  will  have  a 
finite  on-time.  Beam  duration  times  will  be  constrained  as  a  result  of  the 
influence  of  a  number  of  potentially  limiting  factors,  such  as  flow  time,  and 
vibration  to  translation  (V-T)  relaxation  processes  in  the  gas  phase  as  well 
as  on  the  nozzle  v/a  11s. 


54 


For  our  experimental  conditions,  we  anticipate  the  production  of  10 
percent  excitation  of  v  »  1  H2  molecules  in  ?  "slug"  of  gas  with  dimensions 
of  0.8  mm  diameter  and  ~  190  cm  long  moving  at  near  Mach  1  or  1.5  x  10^ 
cm/sec.  It  will,  therefore,  take  1.3  msec  to  flow  this  "slug"  of  qas  out  of 
the  chamber.  This  flow  time  can  be  compared  to  various  potential  loss 
processes  also  expressed  as  characteristic  times,  e.g.,  consider  losses  in  the 
degree  of  excitation  due  to  vibration-to-translation  energy  transfer 
processes,  i.e., 

H2  (v-1)  +  !^H2  +  \\?  +  aE 

where  kUT  at  room  temperature  has  been  determined  to  be  -  1.5  x  10“16 
3  -i 

cm  sec  (see  Figure  21).  This  results  in  a  relaxation  time  of 

VT  °  {^VT  ^2^}  "  ^  msec 


which  is  long  compared  to  the  flow  time,  so  little  de-excitation  due  to 
relaxation  in  the  gas  phase  is  expected. 

Another  consideration  is  relaxation  at  the  walls.  Diffusion  to  the 
walls  can  be  of  two  forms:  mass  transport  and  excitation  transport.  With 
regard  to  mass  diffusion,  a  diffusion  time  is  given  by 


tD 


(0.04)2 

rar~ 


0.01  msec, 


which  is  very  fast  compared  to  the  flow  time  so  the  potential  for 
de-excitation  is  significant.  Alternately,  energy  changing  collisions  may 
provide  a  mechanism  for  even  more  rapid  transport  of  vibrational  excitation  to 
the  walls.  Studies  have  been  done  by  Kovacs,  and  Javan^5^^  for  CC^tOOl) 
diffusion  in  C02  with  the  result  that  diffusion  was  found,  in  fact,  to  bp 
slower  due  to  the  effective  larger  cross  section  involved.  Therefore, 
regardless  of  the  diffusion  mechanism,  H2(v=l)  molecules  will  suffer  tens  of 
collisions  with  the  nozzle  wall  before  exiting.  At  issue  then  is  the  likely 
effect  of  the  wall  for  deactivation.  Studies  by  Heidner  and  Kasper^^  on 
vibrational  relaxation  of  H2(v*l)  also  provided  an  upper  limit  on  the 


58.  Kovacs,  M.A.,  Ph.D.  Thesis,  (Massachusetts  Institute  of  Technology, 
1969). 

59.  Heidner,  R.F.  and  Kasper,  O.V.V.,  Chern.  Phys.  Lett.  15,  179  (1972). 
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Figure  21  Experimentally  Derived  Rate  Constants,  iqo  (M)  as  a  Function  of 
Temperature  T  for  Deactivation  of  H2  by  Various  Collision 
Partners  M 


probability  for  deactivation  of  frivol)  on  collisions  with  pyrex  walls. 

They  found  this  limit  to  correspond  to  a  probability  constant  which 
corresnonds  to  many  thousands  of  encounters.  This  suggests  that  once  formed 
(essentially  instantaneously  by  the  laser  pulse),  the  Hg  metastable 
molecules  may  suffer  numerous  wall  collisions  before  exiting  the  nozzle,  but 
with  little  loss  in  excitation  either  from  gas  phase  or  wall  encounters.  This 
fixes  the  flow  time  to  repetition  frequency  as  the  duty  cycle  or  on-time  of 
the  metastable  beam. 

The  conwercial  Lumonics  laser  which  we  have  used  in  performing  our 
sample  calculations  can  pulse  at  10  pps,  which  allows  a  duty  cycle  of  about 
1  percent,  however,  100  to  1G0O  Hz  devices  have  l-een  built  for  isotope 
separation  applications^6^  and,  therefore,  are  conceptual ly  reasonable  to 
consider.  A  500  Hz  laser  would  give  a  duty  cycle  of  >  50  percent. 

4.  Limiting  Assumptions 

All  of  the  above  discussion  has  considered  the  excitation  laser  sources 
to  be  flat  top  temporal  pulses  of  diffraction  limited  beam  quality.  Clearly, 
tins  is  an  idea  that  is  difficult  to  achieve  experimentally.  This  will  cause 
these  projections  to  be  modified  to  some  extent.  The  fact  that  significantly 
greater  flux  is  projected  than  our  original  objective  gives  us  confidence  that 
such  a  high  flux  is  not  totally  unrealistic.  If  achieved,  it  would  represent 
a  considerable  advance  in  the  current  state  of  the  art. 

We  have  also  defined  our  flux  at  the  exit  port  of  the  nozzle.  Clearly, 
a  "beam''  exiting  a  small  diameter  (0.8  mm)  nozzle  with  a  orivino  pressure  of 
6  torr  will  expand  and  form  a  classic  supersonic  expansion  jet.  Depending  on 
the  application,  one  may  have  to  consider  skimmers  to  redefine  the  beam  at 
some  distance  from  the  exit  nozzle  with  a  subsequent  reduction  in  flux  through 
the  skimmer's  entrance  orifice.  An  entire  literature  has  developed  devoted  to 
understanding  and  characterizing  these  phenomena  and  detailed  calculations  are 
outside  the  scope  of  this  program.  Several  studies  at  similar  pressures  and 
nozzle  diameters  have  been  done  and  the  reader  is  referred  to  these  literature 
sources. 

60.  Pike,  C.,  Avco  Everett  Research  Laboratory,  Inc.  (private  communication). 

61.  See,  for  example,  "Rarefied  Gas  Dynamics,  edited  by  R.  Campargue  (Paris, 

1979)  and  References  therein. 
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Other  potential  limiting  issues  which  we  have  not  fully  evaluated 
include:  gas  breakdown,  plasma  formation  due  to  laser  nozzle  interactions,  and 
clustering  phenomena.  We  do  not  expect  gas  breakdown  to  be  an  issue  in  the 
example  in  that  the  laser  beams  (total  energy  of  80  mO  in  8  nsec  focused 
to  0.8  mm  diameter)  are  focused  to  near  2  GW/cm  .  This  is  below  the  break¬ 
down^2^  of  Hg  and  the  intensity  will  (in  principle)  never  be  achieved  since 
these  beams  will  be  efficiently  converted  to  H2  excitation.  However,  surface 
breakdown  due  to  laser  irradiation  at  solid/gas  interfaces  is  clearly  to  be 
avoided,  since  we  could  expect  plasma  formation  for  these  same  conditions. 

This  requires  a  suitable  alignment  procedure  to  insure  the  laser  light 
propagates  on  axis  centrally  down  the  nozzle.  A  trade  off  of  nozzle  to 
laser  beam  diameter  may  be  necessary  to  avoid  these  difficulties  or  one  may 
consider  going  to  longer  wavelengths  (e.g.,  Ruby,  doubled  yag)  to  reduce  the 
potential  for  plasma  formation. ^ 62 ^ 

Also,  since  the  gas  expands  on  exiting  the  nozzle,  one  may  have  the 
potential  for  hydrogen  cluster  formation,  i.e.,  (H2)n*»  however,  studies^63^ 
to  investigate  the  properties  of  H2  clusters  are  typically  carried  out  at 
source  temperatures  below  80°K  -  not  the  room  temperature  conditions 
considered  here.  Also,  the  influence  of  high  effective  vibrational  tempera¬ 
tures  in  the  hydrogen  beam  (e.g.,  10  percent  excitation  corresponds  to 
Tv  ~  2600*K)  should  further  lessen  the  likelihood  of  this  formation.  In 
fact,  performing  such  experiments  at  low  source  temperatures  with  and  without 
vibrational  excitation  would  be  an  interesting  fundamental  application  of 
this  technique. 


62.  Alcock,  A.J.,  Kato,  K.  and  Richardson,  M.C.,  Optics  Communications  6, 

63.  See,  for  example,  "The  Population  of  1=0  and  L=1  levels,  for  Ortho-Para 
(H2)  in  an  Expanding  Jet"  by  J.  Verberne  and  J.  Reuss  in  Rarefied  Gas 
Dynamics,  edited  by  R.  Campargue  (Paris,  1979)  p.  1049. 
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5.  Other  Applications 

This  technique  for  producing  vibrationaily  excited  molecules  has  found 
limited  application  to  date  in  the  literature.  The  most  significant  has  been 
the  work  of  Oucuing  and  co-workers, '^)  whereby  the  excitation  was  directly 
accomplished  using  a  high  pressure  cell  to  generate  Stokes  radiation 
followed  by  the  "amplifying  cell"  or,  in  more  recent  work, ^^generalizing 
the  technique  by  using  two  dye  lasers  pumping  Ng  and  Og. '  In  these 
experiments,  fractional  populations  were  purposely  kept  low,  so  that  the 
process  of  interest  could  be  clearly  studied,  i.e.,  vibration  to  translation 
energy  transfer. 

This  technique  has  also  round  application  in  chemical  kinetic  studies 
designed  to  evaluate  the  effect  of  vibrational  energy  in  overcoming  chemical 
barriers,  e.g.,  the  reaction 

09  +  H,  DH  +  HD 
(661  i  c 

was  investigated'  '  using  this  technique  to  ascertain  ’f  localized 
vibrational  energy  was  effective  in  overcoming  the  >  25  kcal/mole  activation 
energy. 

It  is  easy  to  envision  other  reactions  with  applications  in  isotope 
separation,  chemiluminscence  studies,  lasers  and  astrophysics. 


64.  Queuing,  J.,  Joffrin,  C.  and  Coffinet,  J.P.,  Optics  Comm.  2,  245  (1970). 

65.  Frey,  R.,  lukasik,  0.  and  Ducuing,  J.,  Chem.  Phys.  Lett.  14,  514  (1972). 

66.  Bauer,  S.H.,  Lederman,  D.M. ,  Resler,  E.L.  and  Fisher,  E.R. ,  Int.  J.  Chem. 
Kinetics  5,  93  (1973). 
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B.  OPTICAL  PUMPING  METHODS  FOR  THE  PRODUCTION  OF  HIGH  DENSITY,  HIGH- 
SELECTIVITY  METASTABLE  SOURCES 
1.  Introduction 

The  basic  concept  is  to  use  several  optical  pumping  sources  ( i . e. ,  fluo¬ 
rescence  sources,  lasers,  or  combinations  of  both)  to  access  metastable  levels 
of  atoms  initially  in  the  ground  state.  Optical  pumping  has  the  advantage  of 
being  a  selective  excitation  method,  as  opposed  to  discharge  or  e-beam  excita¬ 
tion,  in  which  the  free  electrons,  as  well  as  other  transient  species  formed 
during  the  pumping  process,  tend  to  destroy  the  metastables. 

The  metastable  rare  gases  are  of  particular  interest  both  to  basic  kine¬ 
tics  studies  as  well  as  to  the  physics  of  high-power  lasers  (e.g.,  the  rare 

gas-halides).  Unfortunately,  even  the  lowest  of  the  metastable  rare  gas  en- 

3  0 

ergy  levels,  the  Xe*(  P2)  with  an  excitation  energy  of  8.3  eV,  is  out  of 

reach  of  existing  tunable  lasers.  We  have,  therefore,  devised  a  combined 

flashlamp  +  laser  optical  pumping  scheme  to  produce  a  beam  of  metastable  xenon 

1*7  —3 

at  typically  [Xe*]  ~  10  cm  and  with  a  selectivity  (i.e.,  percent  of 
metastables  relative  to  ground-state  atoms)  of  -  10  percent.  It  should  be 
noted  that  such  a  metastable  source  would  be  well  beyond  the  existing  state  of 
the  art. 

Atomic  krypton  at  a  pressure  of  ~  1-2  atm,  pumped  by  a  high-energy 

2 

(-  450  kV),  intense  (~  20  A/cm  )  e-beam,  emits  a  strong  UV  band  centered  at 
X  ~  1460  A  and  with  FWHM  ~  138  A;  this  radiation  is  attributed  to  the  Kr~ 

/n  ro\  ★ 

excimer. v  *  '  The  present  scheme  is  to  use  this  Kr2  source  as  a  "flashlamp" 

to  pump  the  Xe  resonance  line  at  \  »  1470  A.  The  energy  levels  of  atomic  Xe 
and  the  important  kinetic  processes  are  shown  schematically  in  Figure  22. 
Briefly,  the  ground  state  Xe  (5p°  Sn)  atoms  are  excited  by  the  Kr0 

C  *5  A  VJ  L 

flashlamp  to  the  Xe*  (5p  6s  P, )  resonance  level.  These  excited  atoms  are 
subsequently  transferred  with  a  tunable  laser  to  one  of  the  ten  Xe*  (5p  6p) 
states  the  ("2p6"  level  is  preferred)  and  then  spontaneously  decay  to  the 
metastable  Xe*  (5p56s3P2)  state,  which  has  a  natural  lifetime^9^ 

67i  Hoff,  P.W.,  Swingle,  J.C.,  and  Rhodes,  C.K.,  Appl.  Phys.  Lett.  23,  245 
(1973). 

68.  Duzy,  C.  and  Boness,  M.J.W.,  IEEE  OQE  16,  640  (1980);  Boness,  M.J.W.  and 
Duzy,  C.,  "Development  of  Rare  Gas  Photolytic  Drivers,"  Final  Technical 
Report,  Dec.  1979;  private  communication. 

69.  Small-Warren,  N.E.  and  Chiu,  L.Y.C.,  Phys.  Rev.  All,  1777  (197r\ 
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of  ~  15U  sec.  The  various  processes  indicated  in  Figure  22  will  be  discussed 
in  detail  below. 

A  conceptual  view  of  the  apparatus  is  shown  in  Figure  23.  It  consists 
of  a  coaxial  arrangement,  with  an  outer  cell  filled  with  Kr  at  -  1-2  atm  pres¬ 
sure  and  inner  cell,  made  of  MgF?  to  pass  the  1470  A  radiation,  containing 

13  ^  -3 

the  Xe  beam  at  a  density  -  10  °  cm  and  with  a  beam  diemeter  -  1  nm.  Ihe 
laser  radiation  propagates  along  the  axis  of  the  inner  cell.  A  transverse 
e-beam  irradiates  the  Kr  in  the  outer  cell,  forming  the  Kr,,  excimer, 
which  then  optically  pumps  the  Xe  in  the  inner  cell,  as  described  above.  A 
similar  system  (without  the  tunable  laser)  excists  at  AERL,  and  has  yielded  a 
measured^68^  fluence  of  -  100  mJ/cm^  impinging  upon  the  central  tube  with 

C  O 

a  pulsewidth  of  ~  1  nsec,  giving  a  pump  power  of  I  ~  10  W/cm  in  the 
* 

Kr^  band  of  interest.  This  is  clearly  a  very  intense  "f lashlamp,"  and  is 

well-suited  for  our  purposes.  One  caveat  to  be  interjected  at  this  point  is 

that  krypton  gas  usually  contains  trace  amounts  of  xenon  as  an  impurity.  This 

is  clearly  intolerable  for  the  present  scheme,  since  the  xenon  resonance  line 

* 

would  absorb  just  that  portion  of  the  Kr,,  band  that  is  needed  to  pump  the 
coaxial  Xe  beam.  It  is,  therefore,  necessary  to  obtain  pure  (i.e.,  Xe-free) 
krypton.  Should  this  not  be  technically  feasible,  the  present  calculations 
can  be  regarded  as  a  generic  example  of  an  optically-pumped  metastable  source, 
illustrating  the  relevant  principles  and  constraints. 

2.  Results  and  Discussion 
a.  Optical  Pumping  Rates 

We  now  consider  in  more  detail  the  optical  pumping  scheme  shown  in 
Figure  22.  In  order  to  model  the  system  we  must  first  calculate  the  param¬ 
eters  indicated  in  the  figure. 

•k 

We  begin  with  the  Krg  flashlamp  pumping  rate,  Rp.  At  the  den¬ 
sities  of  interest,  the  Xe  resonance  line  is  Doppler  broadened,  with  a  FWHM 
given  by 


2v  [2  In  2  k7  =  0.073  cm-1 

d  =  r*  y  vr 


(10) 
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Ihe  oscillator  strength  of  the  Xe(1SQ  »  3P°)  resonance  transition  is^ 

f  *  0.22,  giving  an  absorption  cross  section  at  line  center  of  an  =  2.5  x 
—12  2  u 

10  cm  .  In  order  for  the  pump  photons  to  be  effectively  absorbed,  we 

must  have  the  optical  depth  q  og[Xe]r  >  1.  Taking  a  beam  radius  of 

r  b  0.5  mm  gives  [Xe]  >  8  x  10  cm-3.  To  obtain  high  metastable  selec- 

13  -3 

tivity  and  reasonably  uniform  deposition,  we  take  [Xe]  =  10  cm  .  It  is 

of  interest  at  this  point  to  compare  the  total  number  of  Xe  atoms  with  the  to- 

5  2 

tal  number  of  pump  photons.  The  pump  flux  is  I  =  10  W/cm  incident  on  the 

pO  p 

beam,  giving  <L  =  I„/hv  *  7.4  x  10“  photons/ sec-cm.  However,  the  total 
P  P  .  3-1 

pump  bandwidth  is  138  A  or  rp  ®  6.5  x  10  cm  ,  while  only  those  photons  within 
approximately  the  Doppler  width  will  be  effectively  absorbed.  Thus,  the  "use¬ 
ful"  pump  flux  is  only 


ft)- 


photons 

2 

sec-cnr 


corresponding  to  I'  =  1.1  W/cm  .  We  are,  therefore,  only  utilizing 
C  P 

~  10  of  the  available  pump  intensity  due  to  bandwidth  limitations.  The 
total  number  of  useful  photons  per  unit  length  of  the  beam  is  then  given  by 


np  =  ^p  Tpulse^nr^  =  2,6  x  10  Photons/cm* 


while  the  total  number  of  ground-state  Xe  atoms  per  unit  beam  length  is 


NXe  =  =  x  aton,s/cm'  (13) 


70.  Ganas,  P.S.  and  Green,  A.E.S.,  Phys.  Rev.  A4,  182  (1971),  and  references 
contained  therein. 
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The  number  of  pump  photons  available  is  therefore  sufficient  to  bleach  the  Xe 
atoms;  however,  the  pump  rate  and  pulsewidth  are  such  that  bleaching  does  not 
occur  (see  below).  The  pump  rate  is  given  by 


/m  4>  i  r  “ 

0  r!  °d*  *  7^  *ref  *  2‘3 


x  106  sec-1 


where  r0  -  2.818  x  10~  cm  is  the  classical  electron  radius  and  all  of  the 
other  quantities  have  been  defined  earlier.  Finally,  knowing  f,  we  can  calcu¬ 
late  Ap  the  spontaneous  decay  rate,  from 


2 

8*  rQc 
e 


2.3  x  108  sec-1 


o  n 

The  next  step  is  to  consider  the  transfer  of  excitation  from  the  Xe*(  Pi) 

3  0  A 

resonance  level  over  to  the  Xe*(  P2)  metastable  state.  As  discussed  briefly 

in  the  Introduction,  this  is  accomplished  via  one  of  the  ten  Xe*(5p56p)  levels 

5 

with  a  tunable  laser.  The  particular  choice  of  5p  6p  level  involves  selec¬ 
tion  rules,  branching  ratios,  and  transition  rates.  The  2p^  state  (see  Ref. 

71  for  notation)  appears  to  be  ideal:  it  has  branching  ratios  of  0.7  to  the 

30  30  (72) 

P2  state,  0.3  to  the  P^  state,  and  no  branching  to  any  other  level.  ' 

The  transition  wavelengths,  A-f actors  (A2  and  Aj  of  Figure  22)  and  absorp¬ 
tion  oscillator  strengths  are  tabulated  below:' 


71.  Moore,  C.E.,  Atomic  Energy  Levels.  U.S.  Nat.  Bur.  Stand.  Circ.  No.  467 
(U.S.  Dept,  of  Commerce,  Washington,  DX.,  1949),  Vol.  III. 

72.  Setser,  D.W. ,  private  communication. 
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lransition 


2pg  -  3P°  8952 

2pg  -  3P°  8232 


AUec”1) 

f 

9.6  x  106 

0.12 

2.2  x  107 

0.13 

Assuming  the  two  transitions  to  be  Doppler-broadened,  the  relevant  cross 

sections  are  oq(3P®  »  2pg)  «  oq^P®  >  2pg)  a  8.3  x  10-12  cm2.  In  order 

to  obtain  high  selectivity,  we  would  like  the  2p,  •>  3P?  transition  to  be 
"  '  b  c 

optically  thin,  so  that  the  metastable  state  acts  as  a  true  excitation  sink. 
Consistent  with  this  picture,  we  have  neglected  reabsorption  of  the  8232  A 
photons  in  our  model  (see  Figure  22).  Ihe  requirement  that  the  8232  A  transi¬ 
tion  be  optically  thin  leads  to  a  constraint  on  the  metastable  density:  tq  b 
°0[Xe*(3P;>)]r  <  1,  giving  [Xe^pjij)]  <  2.4  x  1012  cm-3.  Finally,  the 
rate  for  laser  excitation,  R^,  is  given  by 


,(3P?  *  2  P, 


>1 

K 


3.7  x  10  I,  sec 


-1 


(16) 


where  It  is  the  intensity  in  W/cm2  of  the  tunable  laser  at  8952  A.  The 

L  y 

reverse  rate,  R^,  is  given  by  detailed  balance  to  be  RLQ  b  2.2  x  10 
IL(sec-1). 

The  remaining  processes  involve  photoionization  from  the  excited  states 
due  to  the  short-wavelength  pump  radiation.  It  should  be  note  that  photo¬ 
ionization  involves  a  broadband  continuum  as  the  final  state,  so  that  the  en¬ 
tire  pump  bandwidth  is  effective  in  this  case,  as  opposed  to  the  resonance  line 

_5 

absorption  process  which  utilized  only  ~  10  of  the  pump  bandwidth.  For- 

3  0  3  0 

tunately,  the  photoionization  cross  sections  from  the  Xe*(  Pi)  and  Xe*(  Po) 

(73)  1  3  n  c 

levels  are  very  small  at  the  wavelength  of  interest/  '  with  oDH(  Pi)  ~ 

3  0  P  rn  1 

°PH'  ~  3  x  10  cm,  giving 


73.  Ouzy,  C.  and  Hyman,  H.A.  Phys.  Rev.  A22,  1878  (1980). 
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°PH*p 


2.2  x  103  sec**1 


Ihe  value^73^  for  opH(2pg)  «  10“*8  cm2,  gives  RpH2  ~  7.4  x  10^  sec"1-. 

b.  Modeling  Calculations 

We  have  numerically  modeled  the  system  shown  in  Figure  22  using  the 
rates  derived  above.  In  Figure  24  we  have  plotted  the  metastable  density 
normalized  to  the  initial  ground-state  density  as  a  function  of  1^,  the  in¬ 
tensity  of  the  tunable  laser.  Ihe  fractional  metastable  density  "saturates" 
at  a  value  of  0.26  (i.e.,  bleaching  does  not  occur,  as  indicated  earlier). 
Starting  with  an  initial  Xe  density  of  [Xe]Q  *  10  cm' ,  this  gives  a 
maximum  metastable  density  of  2.6  x  10*2  cm-3.  However,  from  our  earlier 
arguments,  we  required  [Xe*(  Pi>)]  <  2.4  x  10  cm  for  the  transition 
to  be  optically  thin.  Thus  our  model,  neglecting  reabsorption,  is  only  valid 
at  the  lower  fractional  excitation  values.  A  realistic  expectation  would  be  a 

fractional  density  of  [Xe*(  Pn)]/[Xe]n  »  0.1  giving  a  metastable  density  of 

12-3  ^  ^  2 

about  10  cm  .  This  would  require  ~  1  W/cm  of  \  =  8952  A  laser  in¬ 
tensity  over  a  1  psec  pulsewidth.  We  may,  however,  require  somewhat  more  in¬ 
tensity  if  the  laser  cannot  be  narrowed  to  the  Doppler  width  of  the  8952  A 
transition,  which  is  ~  0.01  A. 

•ft 

In  order  to  test  the  sensitivity  of  the  results  to  the  Kr?  pump  in¬ 
tensity,  Ip,  we  have  repeated  the  modeling  calculations  for  the  case  Ip  = 

5  x  10^  W/cm2  (i.e.,  one-half  of  the  original  value  given  in  Section  IV-B-2-a). 
The  results  are  shown  in  Figure  24.  In  this  case  the  fractional  metastability 
"saturates"  at  0.14,  as  compared  to  the  earlier  value  of  0.26.  In  general,  the 
magnitude  of  [Xe*{  ) 3/CXe30  for  comparable  laser  intensity,  1^,  is  seen 

to  be  approximately  proportional  to  the  value  of  the  Kr^,  pump  flux. 

In  order  to  calculate  the  linewidth  for  the  Xe  resonance  transition 
[Eq.  (10)],  we  used  the  Doppler  velocity  <vD>  =-^2k7/M,  with  T  =  300°K  as 
the  source  temperature.  For  a  collimated  atomic  beam,  however,  we  may  have 

<vaxial>  55  <v0>>  but  '"transverse*  «  <V‘  ™s  wu,<1  lead  t0  a 

narrow  line  in  the  case  of  transverse  optical  pumping  with  a  coherent  (i.e., 


laser)  source,  where  the  relevant  parameter  is  <v 


transverse 


>.  For  an  in¬ 


coherent  (i.e.,  flashlamp)  optical  pump,  however,  the  photons  impinge  upon  the 
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beam  at  all  angles  and  therefore  utilize  the  axial  compo/unt  of  the  velocity 
as  well  as  the  transverse  component.  An  accurate  calculation  wouW  require  a 
convolution  of  the  spatial  photon  distribution  for  the  particular  geometry  of 
interest  with  the  velocity  distribution  of  the  atomic  beam.  For  the  flashlamp 
pumping  case  considered  here,  we  do  not  expect  this  effect  to  significantly 
alter  our  results  or  conclusions. 

3.  Conclusions 

We  have  shown  that  the  flashlamp  +  laser  optical  pumping  scheme  under 
consideration  should  be  capable  of  producing  a  metastable  xenon  beam,  with 
metastable  density  -  10  cm  and  selectivity  of  ~  10  percent.  This  rep¬ 
resents  a  large  improvement  over  the  current  state  of  the  art.  The  principal 
constraint  on  extending  the  scheme  to  higher  densities  is  the  broad  bandwidth 
of  the  Krg  fluorescence  pump  relative  to  the  narrow  absorption  bandwidth 
of  the  Xe  atomic  resonance  line.  From  formula  (12),  we  found  that  the  number 

of  useful  pump  photons  was  n'_  «  2.6  x  10^  photons/cm  of  beam  lenqth,  which 

P  13  3 

corresponds  to  a  density  of  3.3  x  10  photons/cm  .  Thus,  even  if  there  were 
no  losses  in  the  system  whatsoever,  so  that  every  pump  photon  yielded  a  meta- 

IT! 

stable  atom,  we  could  still  only  achieve  a  metastable  density  of  -  3  x  10 
atoms/cm3.  Generally,  fluorescence  pump  sources  do  not  have  the  intensity 
per  unit  bandwidth  to  achieve  extremely  high  metastable  densities  ( i . e. , 

-  1016  cm"3). 

As  indicated  in  the  Introduction  the  present  scheme  for  accessing  high- 
lying  metastable  energy  levels,  which  are  typically  out  of  the  range  of  exist¬ 
ing  tunable  lasers,  should  be  generally  applicable  to  a  variety  of  systems. 

For  example,  the  Ar?  excimer  emission  band  overlaps  the  krypton  atomic 
resonance  line,^7^  while  the  Xel*£>£band  overlaps  the  mercury  2537  A 
resonance  transition. In  fact,  at  AERL  we  have  used  a  coaxial,  e-beam 
pumped  device,  of  the  type  shown  in  Figure  23,  to  produce  [Hg*]  >  5  x  1014 
cm"3  optically  pumped  by  the  Xel*  fluorescence  handJ73^  Detailed  model¬ 
ing  for  these,  and  perhaps  other  cases,  could  v  -arried  out  in  a  similar 
fashion  to  the  xenon  case. 


"W.  Wrobel,  W.G.,  Rohr,  H.  and  Steuer,  K.H.,  Appl.  Phys.  Lett.  36,  113 
(1980). 

75.  Mandl,  A.  and  Hyman,  H.A.,  0.  Chem.  Phys.  74,  3167  (1981). 
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Finally,  it  should  be  noted  that  although  the  K r*  "flashlamp"  pro¬ 
duces  a  fluorescence  intensity  -  105  W/cm2,  only  -  1  W/cm2  is  useful 
(i.e.,  within  the  absorption  bandwidth).  It  may  therefore  be  worthwhile  to 
examine  the  possibility  of  utilizing  an  intense,  but  unreversed,  xenon  reso¬ 
nance  lamp  to  accomplish  the  "flashlamp"  pumping  step.  Using  a  coaxial  xenon 
lamp  in  place  of  the  e-beam  pumped  Kr2  source  would  greatly  simplify  the 
apparatus,  and  may  therefore  be  desirable  even  at  the  sacrifice  of  some  met¬ 
astable  density.  Similarly,  intense  resonance  lamps  might  also  be  useful  for 
the  krypton  and  mercury  examples  cited  above.  Ihe  operating  characteristics 
of  these  sources  are  described  in  the  next  section  along  with  a  laser  based 
method  for  the  generation  of  1470*  A  xenon  resonance  radiation  via  a  four  wave 
mixing  technique  in  magnesium  vapor. 


4.  Other  Methods  for  Generating  Vacuum  Ultraviolet  (VUV)  Radiation 

a.  VUV  Generation  by  Four-Wave  Mix  in 

We  wish  to  transfer  a  sizeable  fraction  of  a  ground-state  (*S  )  xenon 
.  3  0  0 
beam  into  the  metastable  P2  level.  We  have  established  that  this  may  be 

done  by  optically  pumping  the  Xe  atoms  with  two  light  sources  of  frequency 
6S045.7  cm-1  (1469.6  A  vacuum  wavelength)  and  11167.3  cm-1  (8954.7  A  wave¬ 
length),  respectively.  We  propose  that  both  sources  be  pulsed,  that  the  IR 


source  be  a  dye  laser,  and  that  the  VUV  be  generated  by  four-wave  mixing  of 
two  dye  lasers. 


A  schematic  of  our  proposed  metastable  generator  is  shown  in  Figure  25. 

A  pulsed  laser  simultaneously  pumps  three  dye  lasers  of  vacuum  wavelengths 
8955,  4310  and  4620  A,  respectively.  (Alternatively  the  pump  can  be  a 
frequency-doubled  and  tripled  Nd-YAG  laser.)  Ihe  visible  lasers  are  mixed  in 
a  cell  containing  **  10  torr  of  Mg  vapor  and  He  buffer  gas;  with  the  wavelength 
of  one  laser  (4310  ft)  tuned  to  a  two-photon  resonance  of  atomic  Mg,  the  wave 
mixing  is  enhanced  and  efficient  production  of  the  VUV  (at  the  sum  frequency 
2WJ  +  \>2)  will  occur  (see  Figure  26).  1470  A  VUV  production  by  this  tech¬ 
nique  was  reported  in  1975  by  Wallace  and  Zdasiuk/76^  These  authors 


^67  Wallace,  S.  and  Zdasiuk,  6.,  Appl.  Phys.  Lett  28,  449  (1976). 
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Figure  25  Generation  of  a  Pulsed  Beam  of  Metastable  Xe  Atoms  Using  Four-Wave 
Mixing  in  Hg  Vapor 
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Figure  26  Energy  Levels  of  Atomic  Mg  Relevant  to  Four-Wave  Mixing  VUV 
Generation.  The  Dye  Laser  Frequencies  are  Labeled  and  v£; 
vjjv  *  2vj  +  vg 


reported  a  conversion  efficiency  (vj  to  vuv)  of  -  0.1  percent  for  peak  dye 
User  powers  of  50  kW  (in  a  ~  7-nsec  pulse)  and  dye  laser  spectral  widths  of 
~  0.1  cm-*.  Because  of  pulse  narrowing  which  is  inherent  in  wave  mixing  of 
this  type,  we  can  expect  to  generate  a  3-nsec-wide,  50  W  peak-power  VUV  pulse, 
corresponding  to  a  VUV  energy  of  ~  0.15  yJ  in  a  spectral  width  of  0.3  cnf*. 
This  can  be  easily  focused  into  an  Xe  beam  to  give  a  fluence  of  10  J/cm  ,  a 
fluence  which  is  sufficient  to  saturate  the  Xe  S^-  P ^  transition  (see  Sec¬ 
tion  IV-B-5).  The  VUV  pulse  is  timed  to  coincide  with  the  center  of  the 
7-nsec  IR  pulse  (8955  A)  in  the  region  of  their  spatial  overlap.  Thus,  the  Xe 

atcms  excited  up  to  the  P,  state  will  tend  to  be  excited  by  the  IR  laser 

A  o 

to  the  2p,  level  if  the  IR  laser  spectral  intensity  (W/cm  GHz)  is  «uf- 
0  3  0 

ficiently  large  to  make  the  P^  -  2pg  transition  rate  Rcb  much  greater 
than  the  rate  A^a  (■  2.3  x  10®  S-1)  at  which  the  3P®  level  decays 
spontaneously  (Figure  27).  The  required  IR  spectral  intensity  is  a  relatively 

p 

modest  80  W/cm  GHz,  which  is  easily  achievable  in  a  pulsed  dye  layer,  (see 

Section  IV-B-5).  Because  of  the  favorable  branching  ratio,  most  of  the  2pK 

3  0  ® 

atoms  decay  spontaneously  to  the  metastable  P £  level  after  the  light  pulses 

end.  Spontaneous  decay  of  the  2pg  level  during  the  light  pulses  is  negli¬ 
gible. 

The  metastable  generator  we  have  described  requires  the  following  hard 
ware:  (a)  a  1  MW  peak-power  pulsed  Ng  laser,  such  as  the  Molectron  UV— 1000; 
(b)  three  Hansch-type  transversely  pumped  dye  lasers,  (c)  two  intracavity  dye 
laser  etalons  (one  in  each  visible  dye  laser)  to  narrow  the  spectrum  to 
~  0.1  cm-1  and,  thereby,  enhance  the  four-wave  mixing,  (d)  ?  heat-pipe- type 
oven  for  atomic  Mg  (the  required  temperature  is  ~  850*C),  (e)  various  focusing 
lenses  and  mirrors,  and  (f)  a  VUV-sensitive  solar-blind  photomultiplier  tube. 
No  elaborate  VUV  spectraneter  is  required,  since  the  proposed  wavelength  tun¬ 
ing  may  be  established  by  observing  the  fluorescence  produced  when  the  VUV 
beam  is  incident  on  a  cell  containing  Xe. 
b.  Low  and  High-Pressure  VUV  Lamps 

There  are  two  types  of  incoherent  VUV  sources  that  are  currently  being 
considered:  a  rare  gas  atomic  line  source  and  an  excimer  emission  band 
source.  The  atomic  line  source  is  a  low  pressure  (^  1  torr),  cw,  microwave 
discharge  lamp  that  contains  a  rare  gas  that  generates  the  appropriate  reso¬ 
nance  fluorescence,  i.e.,  a  Xe  lamp  to  pump  the  Xe  resonance  line,  etc.  A 
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Figure  27  Energy  levels  of  Atomic  Xe  Relevant  To  Metastable  (state  m) 

Generation.  The  level  Heights  (about  0)  are  not  drawn  to  scale. 

The  R's  represent  Light  Induced  Transition  Rates,  while  the  A's  are 
Spontaneous  Decay  Rates.  It  is  seen  that  the  Ratio  Rcb^ba 
determines  what  fraction  of  the  VUV-Induced  Population  in  level  b 
is  transferred  to  level  c.  Because  of  the  favorable  branching 
ratio  (Acm  «  2Acb),  most  of  the  population  excited  to  level  c 
decays  to  the  metastable  level  m 
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number  of  these  devices  have  been  built  and  their  characteristics  reported  in 

the  literature. Photon  fluxes  in  the  range  1016-1017  photons/sec 

(2  x  10“3  -  2  x  10"2  W/cm2  for  the  Xe  resonance  radiation)  were  obtained 

using  a  25  mm  diameter  liF  window.  The  relevant  VUV  fluence  can  be  obtained 

from  the  product  of  the  lamp  flux  and  the  transit  time,  tp,  of  the  atom  beam 

across  the  lamp  window.  For  thermal  Xe  beams,  t  *  115  ysec  for  a  2.5  cm 

P  2 

diameter  window.  Using  a  mid-range  value  for  the  flux,  0.01  W/cm  ,  yields  a 
VUV  fluence  ~  1.2  yJ/cm2.  Assuming  that  the  lamp  linewidth  is  dominated  by 
Doppler  broadening  at  low  pressure  (~  0.1  cm"1)  yields  a  spectral  fluence 
-  0.4  yJ/cm2-GHz,  a  value  comparable  to  that  required  to  saturate  the  reso¬ 
nance  line  transition  (see  Section  IV-B-5).  We  recognize  that  the  detailed 
relative  spectral  power  distribution  of  the  lamp  (including  self-reversal) 
will  have  to  be  measured  to  validate  this  linewidth  assumption.  This  type  of 
measurement  requires  a  very  high  resolution  VUV  monochromator  (ax/x  ~  1  x  10"^) 
that  is  not  available  in  most  laboratories.  Further  investigations  will  be 
required  to  determine  if  this  information  is  available  in  the  literature. 

Broadband  Ar*  (1300-1600  A)  and  Kr?  (1500-1800  A)  excimer  emission 
flashlamps  provide  alternative  continuum  VUV  sources  for  pumping  Kr  and  Xe 
resonance  lines,  respectively.  These  sources  utilize  cold  cathode,  pulsed 
discharges  (Tpuise  =  1-2  usee)  at  several  hundred  torr  noble  gas  pressure. 

This  type  of  excimer  flashlamp  is  more  compatible  with  the  types  of  experi¬ 
ments  described  here  than  is  the  e-beam  pumped  flashlamp  described  previously. 
Although  a  great  deal  is  known  about  the  relative  spectral  bandwidth  of  these 
sources  (ax  -  150  A  FWHM) ,  little  information  is  available  concerning  absolute 
spectral  power  output  of  these  lamps  n  comparison  with  the  absolute  fluores¬ 
cence  efficiencies  available  for  e-beam  pumped  flashlamps.  We  estimate  that 
equivalent  blackbody  temperatures  greater  than  20,000*K  would  be  required  to 
generate  the  necessary  spectral  fluence  to  saturate  the  Doppler  broadened  res¬ 
onance  line  transition. 

The  requirements  on  the  infrared  laser  pump  depend  to  some  extent  on  the 
type  of  incoherent  VUV  excitation  source  that  will  be  employed.  A  cw  tuneable 
dye  laser  would  be  required  if  the  low  pressure  atomic  line  source  is  used. 


Stief,  L.J. 


>1.  Optics  4,  1674  (1965). 
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On  the  other  hand,  the  pulsed  VUV  flasnlamps  would  require  flashlamp  pumped  IR 
dye  lasers  that  have  about  the  same  pulsewdith  as  the  VUV  flashlamp.  Tuneable 
radiation  in  the  7000-8000  A  region  is  required  to  pump  the  s-p  transitions, 
b»c,  in  figure  27  with  a  specti  1  flux  of  ~  80  W/cm  -6Hz.  There  are  at  least 
two  reports  of  cw  tuneable  infrared  radiation  from  krypton  ion  pumped  dye 
lasers  using  carbocyanine  and  oxazine  dyes.^78,79)  Tuned  output  powers  in 
the  range  0. 1-1.0  W  were  obtained  with  a  1  A  linewidths.  The  required  spec- 

O 

tral  flux  (-  100  W/cm  )  can  then  be  obtained  by  focusing  to  a  1  mm  diameter 
spot  (see  Section  IV-B-5).  Similar  results  can  be  expected  from  flashlamp 
pumped  dye  lasers. 

5.  Calculation  of  Required  VUV  Fluence  and  IR  Flux 

We  wish  to  calculate  the  following  two  quantities:  (see  Figure  27) 

(i)  the  VUV  spectral  fluence  required  to  saturate  the  resonance  transition, 

Sq  »  Pj  (A*b)  in  the  rare  gases  and  (ii)  the  IR  spectral  intensity  re¬ 
quired  to  transfer  most  of  population  to  the  higher  p  levels  (b»c). 

Both  requirements  may  be  estimated  by  noting  that  the  transition  rate  R  for  a 
transition  of  wavelength  A  and  oscillator  strength  f  induced  by  a  spectral  in- 

p 

tensity  dl/dv  (erg/cm  Hz)  on  resonance  is 


R 


dl 

ft 


A  transition  is  saturated  if  R  >1,  where  is  the  laser  pulsewidth.  Requir- 

_5 

ing  Rba  t  »  1  for  the  a>b  transition  (xba  =  1.5  x  10  cm,  fba  *  0.22) 
implies  a  minimum  VUV  spectral  fluence  dF/dv  of 

7T  =  ■— T  V--  -  2.3  x  ID"9  erg/cm2  Hz 
37  *reAbafba 


*  0.23  uJ/cm2  GHz 


7 W.  Romanok,  K.M.,  Hildebrand,  0.  and  Gobel,  E.,  Opt.  Comm.  21,  16  (1977). 

79.  Donzel,  A.  and  Weisbuch,  C.,  Opt.  Comm.  J7,  153  (1976). 
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Similarly,  we  calculate  the  minimum  spectral  intensity  of  the  IR  laser 
from  the  requirement  that  Rcb  £  Aba.  7o  be  specific,  we  require  that 
Rcb  ■  5Aba.  With 

xcb  -  9.0  x  10"5  cm,  fcb  -  0.12,  Aba  -  2.3  x  108  s"1, 


we  obtain 


5h 


*rexcbfcb 


ha 


0.81  erg/cm2/s  Hz 


81  W/cm2  GHz 


C.  E-BEAM-SUS7AINED  DISCHARGE  EXCI TA7I0N  OF  METASTABLE  N2(A32*) 

7he  objective  of  this  phase  of  the  program  is  to  assess  the  feasibility 
of  utilizing  an  e-beam  sustained  discharge  as  a  high  flux  molecular  beam 
source  of  nu-tastable  N2(A  eu)  states,  hereafter  abbreviated  A  or  N2(A) 
states.  The  emphasis  in  this  part  of  the  study  is  focused  primarily  on  the 
energy  transfer  and  chemical  kinetic  modeling  aspects  of  high  flux  A-state 
production  using  a  computer  kinetics  code.  The  time  dependence  of  N2(A) 
states  as  well  as  a  number  of  other  species  is  investigated  as  a  function  of 
e-beam  current  density  and  electric  field  strength.  We  also  briefly  discuss  a 
method  for  extracting  A-states  from  the  discharge. 

The  basic  idea  is  to  produce  A-states  by  e-beam  sustained  discharge 
pumping  a  high-speed  flow  (flow  speed  ~  10^  cm/sec)  of  pure  N2  at  near 
atmosphere  pressure  and  then  forming  a  high  flux  molecular  beam  source  by  ex¬ 
panding  a  small  portion  of  the  'low  through  a  supersonic  nozzle/skimmer  sam¬ 
pling  orfice  into  a  high  vacuum  chamber.  Using  this  technique  we  believe  that 

iq  _p  _i 

an  A  state  flux  of  at  least  10  cm  sec  can  be  obtained  at  the  throat 
of  a  sonic  orifice  for  modest  energy  loadings  (30  J/l-atm).  Our  modeling 
calculations  also  indicate  much  higher  flux  levels  can  be  achieved  at  higher 
energy  loadings  (500-1000  J/l-atm).  However,  these  results  require  more 
detailed  kinetic  information  regarding  translational  heating  rates  and 
temperature  dependences  of  a  number  of  important  energy  transfer  processes 

before  an  adequate  assessment  can  be  made. 

There  are  two  distinct  advantages  which  favor  the  use  of  the  e-beam  sus¬ 
tained  discharge  over  the  self-sustained  discharge:  (1)  uniform,  volumetric 
excitation  at  high  gas  densities,  and  (2)  independent  control  over  both  the 
secondary  electron  density  and  electron  energy  distribution.  The  ionization 
rate  is  controlled  by  varying  the  primary  e-beam  current  while  the  electron 
energy  is  controlled  by  varying  the  sustainer  electric  field  strength.  In 
this  way,  the  electron  energy  distribution  can  be  tailored  to  efficiently  pump 
the  electronic  states  in  N2  while,  at  the  same  time,  maintaining  a  rela¬ 
tively  low  ionization  rate  and,  therefore,  a  stable  discharge  at  high  pressure. 

Or  modeling  calculations  of  A-state  excitation  indicate  that  the  total 
lifetime  of  N2(A)  within  the  discharge  zone  is  of  the  order  of  1  msec. 
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Therefore,  the  high-speed  flow  and  supersonic  nozzle  sampler  will  be  required 
to  extract  the  A-states  from  the  discharge  zone  and  transport  them  to  some 
other  region  where  their  properties  can  be  investigated. 

1.  Overview  of  Modeling  Calculation 

The  time  dependence  of  Ng{A)  state  formation  and  loss  is  obtained  in  a 
two  step  calculation.  First  the  AERL  Boltzmann  computer  code  is  used  to  ob¬ 
tain  steady-state  electron  energy  distributions  f((5jS'  and  average  electron 
excitation  rate  constants,  k^(E/N),  at  a  given  value  of  E/N,  the  pressure 
reduced  electric  field  strength.  These  quantities  are  related  by 

dU  f(U)  *  o  .(U)  U  (18) 

o 

where  <^(U)  is  the  energy  dependent  cross  section  for  the  ith  electron-heavy 
particle  inelastic  process: 

e(U)  +  N2(k)  >  N2U)  +  e  (U  -  aUk>Jl)  (19) 

and  where  the  indices  k  and  label  two  N2  quantum  states  and  aU  is  the 

energy  difference  between  those  two  states.  The  values  of  k.  for  electronic 

excitation  and  ionization  are  then  fed  into  a  kinetics  code  calculation  which 

integrates  the  coupled  set  of  rate  equations  that  describe  all  the  formation 

and  loss  processes  included  in  this  model  for  six  species:  N.(A),  N9(B), 

4  2  c  c 

N2(C),  N(  S),  N(  P)  and  electrons  that  are  formed  in  the  discharge.  The 

output  of  the  kinetics  computer  code  is  then  displayed  in  the  form  of  computer 

generated  plots  of  the  temporal  dependence  of  the  number  density  of  electrons, 

N2(A),  N(4S)  and  N(2P)  states. 

2.  Boltzmann  Code  Calculation 

By  utilizing  the  two-step  calculation  described  above,  we  have  made  the 
implicit  assumption  that  the  time  evolution  of  the  secondary  electron  energy 
distribution  can  be  completely  decoupled  from  the  time  evolution  of  the  number 
densities  of  the  electrons  and  heavy  particles.  This  approximation  is  valid 
in  the  present  set  of  calculations  since  the  time  scale  required  for  the 
secondary  electron  distribution  to  settle  down  to  a  steady-state  value  is  much 
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faster  than  the  time  scale  for  variations  in  the  species  number  densities  in 
the  discharge.  Given  that  the  cross  section  for  some  of  the  important  elec¬ 
tron  energy  loss  processes.  a>"C  in  the  range  10~16  -  1(T*7  cm2  and  that 
the  initial  ground  state  N2  density  is  assumed  to  be  1  amagat  (2.45  x  1019  cm"3), 
we  estimate  that  the  characteristic  time  for  electron  energy  relaxation  is  <  1 
nsec  for  1  eV  electrons,  whereas  it  will  be  shown  below  that  number  density 
variations  always  occur  on  much  slower  time  scales. 

7able  2  lists  the  electron  energy  loss  processes  included  in  the  AERL 
Boltzmann  Code  calculation.^30"®7^  Cross-section  values  from  threshold  to 
100  eV  electron  energy  were  used  to  describe  each  process.  Experimental 
energy  dependent  cross-section  data  were  available  for  each  of  these  processes 
except  for  the  three  metastable  ionization  processes  from  N2(Aeu), 
and  N2(C3vu)  states  and  rotational  excitation  in  ground  state  N,,.  Cross  sections 
for  ionization  from  the  A  state  were  taken  from  the  binary  encounter  calculations  of 
lon-lhat  and  Flannery. ^87^  Ionization  cross  sections  from  the  B  and  C  states^88) 
were  then  scaled  from  the  A  state  results  using  the  relationships: 

•8f(U)  -  C°Ai(U  "  <W>‘3  'V:B)2  (20a) 

°Ci (0)  *  C°Ai(U  -  <ic/1a>i>']  (ia/ic>2  (20b> 

whe-'-e  I  and  are  the  ionization  energy  and  cross  section  for  each 
state  X  (=  A,  8,  and  C). 

Using  this  set  of  cross-section  data,  Boness^89^  has  shown  that  the 
calculated  values  of  N2  drift  velocity,  characteristic  energy  and  first 
Towns!:  i  ionization  coefficients  derived  from  the  AERL  Boltzmann  code  are  in 
good  agreement  with  experimental  measurements  of  these  quantities  for  E/N 
values  up  to  ~  6  x  10-15  V-cm2. 

Figure  28  displays  the  steady-state  electron  energy  distribution  at  two 
different  F/N  values:  4.1  x  10“*8  V-cm2  and  1.02  x  10“*8  V-cm2.  The 
average  electron  energy  <U>  for  each  of  these  distributions  is  1.09  and 
1.77  eV,  respectively,  where 

87!  Ton-lhat,  0.  and  Flannery,  M.R.,  Phys.  Rev.  A15,  517  (1977). 

88.  Hyman,  H. ,  Avco  Everett  Research  Laboratory,  Inc,  private  communication. 

89.  Boness,  J.,  Avco  Everett  Research  Laboratory,  Inc.  private  communication. 
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TABLE  2 

ENERGY  LOSS  PROCESSES  INCLUDED  IN  AERL  BOLTZMANN  CODE 

AU  (eV) 


Refs. 


3.938  x  10"5 
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o  Momentum  Transfer 

¥e's 

o  Rotational  Excitation 

es+N2(j) - ►  e's  +  N2(J') 

o  Vibrational  Excitation 

e$+N2(V=0) - ►  e's+N2(V') 

AV  =  1 

AV  =  2 

AV  =  3 
AV  =  4 
AV  =  5 
AV  =  6 
A  V  =  7 
AV  =  8 

d  Electronic  Excitation 

VV*1  rg) - ►  e"s*Nj(A3tJ) 

- *  e'  s+Nz(B3,rg) 

”  - *  e*  s+N2(W3a„) 

"  - »  e'  s+Ng(B' 3l“) 

"  - *  e‘  S+Nz(a- »g) 

"  - *  e1  5+*2(w\) 

“ - »  e' s4N2(C3»u) 

“  - *e's+N2(EV) 


0.02  81 

0.29  82 

0.58 
0.88 
1.17 
1.46 
1.75 
2.Q5 

2.34  " 


6.17  83 

7.35 

7.36 
8.165 

8.55  « 

8.895  ■' 

11.03 

11.88 
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TABLE  2  (Continued) 


PROCESS 


Dissociation 

es+¥X)“^ 

2N(4S) 

Ionization 

vvx>— - 

N2  +  2es 

es«Hj(A>— - 

Nj  (A2ttu)  +  2e, 

VN2<B)— - 

N+  (A2ttu)  +  2e 

VN2<C)— ' 

Ng  (A2*u)  +  2e; 

4U  (eV) 

Refs. 

9.76 

84,  85 

15.58 

86 

10.79 

87 

9.62 

this  work 

5.99 

this  work 

80.  Englehardt,  A.M. ,  Phelps,  A.V.  and  Risk,  C.6.,  Phys.  Rev.  135, 
Al 566  (1964). 

81.  Chandra,  N.  and  Burke,  P.G.,  0.  Phys.  B  6,  2355  (1973). 

82.  Chandra,  N.  and  Temkin,  A.,  Phys.  Rev.  A  J13,  188  (1976)  and 
references  cited  therein. 

83.  Cartwright,  D.C.,  Traji.iar,  /. ,  Chut ji an ,  A.  and  Williams,  W., 
Phys.  Rev.  A  16,  1041  (1977). 

84.  Winters,  H.F.,  Horne,  O.E.  and  Donaldson,  E.E.,  0.  Chem.  Phys, 
41,  2766  (1964). 

85.  Niehaus,  A.  and  Naturforshung,  Z, ,  7?as  690  (1967). 

86.  Rapp,  0.,  Briglia,  0.0,  and  Englander-Golden,  P.,  0.  Chem. 
Phys.  42,  4081  (1965). 

87.  Ton-That,  D.  and  Flannery,  M.R. ,  Phys.  Rev.  A  15,  517 
(1977). 
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Figure  28  Dotted  Curves  are  Electron  Energy  Distributions  as  Calculated  by 
the  Boltzmann  Code:  E/N  -  4.1  x  10“16  v-cmz  (curve  A);  E/N  * 
1.02  x  10-15  v-cm2  (curve  b).  Solid  lines  are  Maxwell 
Boltzmann  Distributions  for  two  average  energies,  Te  ■  1.09  eV 


ft 


<U>  -  f dU  f(U)  U3/2  (21) 

o 

For  comparison,  the  Maxwell-Bolt2mann  (M-B)  distributions  which  yield  the  same 
average  energies  are  also  displayed.  It  can  be  seen  that  the  M-B  distribution 
is  a  better  approximation  at  the  higher  E/N  value,  but  it  always  tends  to 
overestimate  the  high  energy  tail  and  underestimate  the  low  energy  peak  at 
both  E/N  values.  These  discrepancies  should  be  kept  in  mind  since,  later,  we 
will  estimate  rate  constants  for  superelastic  relaxation  of  the  N2(A)  and 
N(2P)  states  from  the  up-excitation  rate  constants  using  detailed  balancing, 
a  procedure  that  is  only  strictly  correct  for  M-B  distributions. 

Fractional  energy  deposition  into  a  particular  inelastic  channel  is 
obtained  by  integrating  Eq.  (18)  for  each  inelastic  process  and  then  evaluat¬ 
ing  the  expression 

t,  -  alLk./Ek  aU 

1  i  i  z  z  i 

wnere  *  is  an  index  running  over  all  the  inelastic  processes.  The  results  for 

total  vibrational  excitation,  $v;  electronic  excitation  into  A,  B  and 

C  states  dissociation  and  ground  state  ionization  are  displayed  in 

Figure  29  as  a  function  of  E/N.  It  can  be  seen  that  vibrational  excitation 

-16  2 

dominates  for  E/N  £  4  x  10  V-cm  where  the  average  electron  energy  is 
about  1  eV  or  less.  This  is  the  well-known  region  where  the  C02  laser  oper¬ 
ates.  At  higher  E/N  the  vibrational  contribution  decreases  and  electronic 
excitation  becomes  important  once  the  high  energy  tail  of  the  electron  energy 
distribution  surpasses  the  first  electronic  threshold.  This  defines  the  lower 
limit  of  the  E/N  range  that  will  lead  to  significant  production  of  N2(A) 
states.  However,  there  is  also  an  upper  limit  to  the  useful  E/N  range  due  to 
the  proportionately  higher  production  rate  of  N(4S)  atoms,  produced  through 
electron  impact  dissociation  of  Ng.  It  will  be  shown  below  that  N(4S) 
atoms  are  the  dominant  quenchers  of  N2(A).  Therefore,  it  is  important  to 
produce  as  few  N{  S)  atoms  as  possible  while  still  producing  as  many  N2(A) 
states  as  possible.  This  effect  can  be  seen  more  clearly  in  Figure  30  which 
displays  the  ratio,  R,  of  the  sum  of  the  A,  B  and  C  electronic  excitation  rate 
constants  k^.  to  the  dissociation  rate  constant,  kQ 
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E/N  (V-em2)  x  I0{5 


Figure  29  Fractional  Energy  Deposition  in  N2  at  1  atm  and  298*K  vs  E/N 


1 


R  «  Jjk£1/kD  i  *  A,  8  and  C  \22) 

versus  E/N.  Ihc  ratio  approaches  one  as  E/N  approaches  4  x  10"^  V-cm^ 
implying  that  at  least  one  N{4S)  atom  is  available  to  quench  each  N2(A). 

On  the  other  hand»  R-values  in  the  range  of  10-100  are  available  at  somewhat 
lower  E/N  where  significant  N/>(A)  production  may  still  be  realized.  As  long 
as  the  N2(A)  states  can  be  extracted  from  the  discharge  before  each  N(^S) 
can  be  cycled  through  10-100  quenching  events,  much  of  the  initial  N2(A) 
state  population  may  still  be  extracted  from  the  discharge  as  will  be  shown  in 
the  next  section. 

3.  Kinetic  Model 

Ihe  following  kinetic  scheme  was  chosen  to  inode!  the  production. 
Ionization  and  Recombination  Processes 


(1)  ep  +  N2 — ►N2+  +  ep  +  e$ 


(2)  e$  +  Ng  (X) 


■N  +  K  +  e  +  e$ 
— -  N2+  *  2es 


(3)  e  +  * 


products 


A-  B-  and  C-State  Excitation 
{4)  e$+  N2(X) - -N2(A}  +  es 

(5)  e$  +  N2(X)—  N2(8)  +  es 

(6)  e$  +  N2(X)  — -  N2(C)  +  es 

(7)  N2(B)  +  N2(X)  — ►  N2(A)  +  N2(X) 

(8)  N2(C)  +  N2(X)  — -N2(B)  +  N2(X) 


(9)  N2(B)  — ►  N2(A)  +  (hvj  -  10,600  A) 

(10)  N2(C) — *-  N2(B)  +  (hv2  .  3371  A) 


N«  Dissociation 


(11)  es  +  N2 — *-N(4S)+N(4S) 

A-  B-  and  C-State  Loss  Processes 

(12)  2  N2(A)— *-N2(B)  +  N2(X) 

(13)  2  N2(A) — ^N2(C)  +  N2(X) 

(14)  N2(A)  +  es — ►  N2  (A)  +  2es 

(15)  N2(B)  +  e$  — N2+(A)  +  2es 

(16)  N2(C)  +  e$  — -N2+(A)  +  2es 

(17)  N2(A)  +  N(4S)  — *-N2(X)  +  N(2P) 


(20)  N(2P)  +  e„ — ►  N(4S)  +  ec 


n(4s) 


4  {+  N2> 
N(4S) - V 


n2*  +  n2 


No(X)  +  No  (X) 


where  N2(X)  is  the  N 2  ground  electronic  state. 

The  N2(A)  states  are  produced  via  the  direct  secondary  electron  impact 
process  (4)  and  through  cascade  from  the  higher  lying  N2(B)  and  N2(C)  states, 
processes  (7),  (8),  (9),  and  (10).  The  primary  N„(A)  loss  process  is  due  to 

4  c.  2 

the  electronic-electronic  energy  transfer  to  N(  S)  atoms  to  produce  N(  P) 

p 

states,  process  (17),  The  N(  P)  states  can  then  relax  via  processes  (18), 
(19),  ana  (20)  allowing  recycling  of  N(4S)  population.  N(4S)  quenchers 
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are  produced  through  both  primary  and  secondary  electron  impact  dissociation 
processes  (1)  and  (11)  and  are  permanently  lost  only  through  recombination 
(21).  Other  molecular  N2  triplet  manifold  loss  processes  are  metastable 
ionization  (14),  (15),  and  (16)  and  superelastic  loss  (17a).  We  now  describe 
how  the  values  for  the  source  functions  and  rate  constants  were  chosen  for 
each  process. 

The  ionization  rate  by  primary  electrons  was  obtained  from  the  energy 
deposition  tables  of  Berger  and  Seltzer^  and  the  energy  partitioning  re¬ 
sults  available  from  the  radiation  chemistry  literature The  total 
e-beam  deposition  rate  in  N2  is  given  by  S  (W  cm-3) 


(* 


PN, 


Jeb 


(23) 


where  Jeb  is  the  primary  e-beam  current  (in  A/cm2),  is  the  density 
of  N2  gas  and  (1/p  3E/3X)^  is  the  stopping  power  at  a  particular  primary 


electron  energy.  For  several  hundred  keV  primary  electrons,  Ref.  90  gives 
(1/p  3E/3X)n  *  2  MeV  -  cm^/gr  yielding  a  total  energy  deposition  rate  at 
1  atm  and  298*K 


2400  Jeb  W/cmw 


(23a) 


The  ionization  rate,  S.  (cm-3  sec"1)  is  then  given  by: 


Sj  -  (S/Wj)  fj  (24) 

where  Wi  is  the  ion  pair  formation  energy  (eV  per  ion  pair  formed)  for  both 
N  and  N2  formation  and  f..  is  the  fraction  of  e-beam  energy  converted 
into  ions.  Reference  91  indicates  W.  =  34.5  eV  and  f.  ^  0.5  yieldinq  S. 

20-3  1  i  i 

=  2.28  x  10  0eb  ions/cm  -sec.  The  secondary  ionization  rate  constant, 
k2,  is  obtained  from  the  Boltzmann  code  and  is  displayed  in  Table  3. 

The  dominant  electron  loss  processes  in  pure  N2  are  assumed  to  be  vol¬ 
umetric  electron  ion  recombination  processes  involving  and  clusters 
formed  in  the  association  reactions: 

"9^*  Berger,  M.J.  and  Seltzer,  S.M*,  Tables  of  Energy  Losses  and  Ranqes  of 
Electrons  and  Positrons.  NASA  SpecTaT  Pub  1 .  3Ui2,  1964. - - 

91.  Willis,  C.  and  Boyd,  A.W.,  Int.  J.  Radiat.  Phys.  Chem.  8,  71  (1976). 


TABLE  3 

EXCITATION  RATE  CONSTANTS  VS  E/N 


RATE  CONSTANT* 

E/N  (X  1015  V-cm2) 

0.41 

0.61 

1.02 

k2 

1.24"20 

9.92“17 

1.47"13 

k4 

2.0~14 

1.61"12 

5.12”11 

k5 

1.68"14 

1.74"12 

6.66"11 

k6 

2.53“16 

1.04~13 

1.24"11 

kll 

5.66"16 

1.5l“13 

1.36"11 

k14 

2.56"17 

2.15~14 

4.95~12 

k15 

1.97-15 

4.29“13 

3.22"11 

k16 

3.68"13 

2.85"11 

8.98"10 

a~N  »  a  x  10"N  cm3/ sec 


*NOTE:  refers  to  rate  constant  for  ith  process  in  the  kinetic  model. 
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(25) 


N+  +  N2  +  N2  — N3+  +  N2 
N2+  +  N2  +  N2  — ►  N4+  +  N2  (26) 


At  these  pressures  the  time  scale  for  clustering  is  orders  of  magnitude  faster 
than  the  expected  recombination  lifetime  (<  10"6  sec)  justifying  the  assump- 
tion  that  N^  and  N^  are  the  dominant  positive  ions  i->  an  e-beam  sus¬ 


tained  discharge.  Ihe  global  recombination  rate  constant  in  pure  N2  has 

been  measured  by  D.  Oouglas-Hamilton^9^  over  the  t/N  range:  10“17  <  E/N  < 

16  2  —  — 
2.2  x  10  V/cm  .  No  attempt  was  made  to  deconvolute  the  separate  recom¬ 
bination  coefficients  for  N^  and  in  these  measurements.  We  fit 


these  measurements  to  the  following  exponential  functional  forms: 


for  1  x  1U-17  <  E/N  <  4  x  10~17  V-cm2: 

k 3 ( E / iN )  -  exp  (-13.28  -  (7.32  x  1016  E/N)) 

for  4  x  Kf17  <  E/N  <  2.2  x  10"16  V-cm2: 

k3(E/N)  «  exp  (-16.15  -  (1.32  x  1015  E/N)) 

and  the  latter  expression  was  then  linearly  extrapolated  (on  semilog  paper)  to 
the  E/N  range  of  interest  here:  0.4  -  1  x  10-*5  V-cm2. 

Ihe  electron  pumping  rate  constants  into  the  A,  B  and  C  states  (k^, 
k5,  kg)  are  listed  in  Table  3  at  three  E/N  values.  Ihe  dominance  of  A  and 
B  state  excitation  is  in  accord  with  other  Boltzmann  cods  results  in  pure 
Ng.'93)  The  quenching  of  N2(C)  states  by  N2  dominate  the  cascade  from 
C  to  B  (kg  =  1.5  x  10"11  cnr/sec  and  k1Q  =  2.2  x  107  sec-*)^94^ 

Similarly,  cascade  from  B  to  A  is  also  dominated  by  col li sional  quenching 
(k?  -  2.0  x  10‘12  cm3 /sec  and  kg  -  1.1  x  105  sec_1).(94) 


92.  Oouglas-Hamilton,  O.H.,  J.  Chem.  Phys.  58,  4820  (1973). 

93.  Judd,  0.,  IEEE  Journal  of  Quantum  Electronics  QE-12  78,  1976. 

94.  Hill,  R.M.,  Gutcheck,  R.A.,  Huestis,  D.L.,  Mukherjee,  D.  and  Lorents, 
D.C.,  "Studies  of  E-Beam  Pumped  Molecular  Lasers,"  SRI  Tech.  Rept. 
MP— 74— 39,  Tech  Rept.  No.  3  Contract  No.  N00014-72-C-0478,  Stanford 
Research  Institute,  Menlo  Park,  California  94025. 
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There  are  two  contributions  to  the  dissociation  yield  of  fu  S): 
secondary  and  primary  electron  impact.  There  are  two  measurements  of  the 
secondary  electron  impact  dissociation  yield  of  Ng.  Winters^)  has 
inferred  the  total  cross  section  for  process  (11),  the  dissociative  excitation 
process  (27a),  and  dissociative  ionization  process  (27b) 

e. 


*  N2_ 

N*  +  N 

(27a) 

+  n2~ 

-*-N+  +  N* 

(27b) 

(where  N*  is  an  excited  nitrogen  atom)  using  a  complex  surface  absorption  ex¬ 
periment  that  required  making  pumping  speed  measurements  and  assuming  values 
for  the  probability  that  N,  N  and  N*  would  stick  to  a  collection  surface. 
These  results  are  in  surprisingly  good  agreement  with  the  (more  straight¬ 
forward)  molecular  beam  measurements  of  Niehaus.^5^  These  latter  measure¬ 
ments  detect  all  N  atoms  but  suppress  all  N+  ions.  Zipf^^^  has  compared 
the  data  of  both  of  these  studies  and  his  results  suggest  that  (11)  and  (27a) 
dominate  (27b).  The  rate  constants  derived  from  these  cross-section  data  are 
displayed  in  Table  3.  The  source  function  for  primary  electron  impact  disso¬ 
ciation  was  obtained  from  the  Berger  and  Seltzer  tables^)  and  the  N0 

(91)  2 

radiation  chemistry  literature'  '  in  a  manner  analogous  to  ionization 

source  function  described  above.  The  energy  required  to  form  a  N(4S)  atom, 

W  ,  is  32.9  eV  and  f  ~  0.5  yielding  an  atom  source  function  rate,  S  , 

a  20a  -3  -1  a 

equal  to  2.28  x  10  atoms  cm  sec  . 

There  is  still  much  disagreement  as  to  the  rate  constant  for  the  pooling 

reactions  (12)  and  (13).  Hays  and  Oskam^6^  have  measured  a  rate  constant 

-9  3 

for  total  loss  out  of  the  A  state,  1.4  x  10  cm  /sec,  and  a  specific 


C  state  formation  rate  constant 


2.6  x  10“^  cm^/sec. 


Stedman  an^ 


Setserv^  and  Zipf^®)  have  measured  2  x  10“^  cm^/sec  for  a  C  state 

production  rate  constant  and  Zipf^9®^  has  obtained  a  total  loss  rate  con- 
-10  3 

stant  =  3  x  10  cm  /sec.  The  Stanford  Research  Institute  kinetics 


95.  Zipf,  E.C.  and  McLaughlin,  R.W.,  Planet.  Space  Sci.  £6,  449  (1978). 

96.  Hays,  G.N.  and  Oskam,  H.J.,  0.  Chem.  Rhys.  j>9,  6088  (1973). 

97.  Stedman,  D.H.  and  Setser,  D.W.,  J.  Chem.  Phys.  50,  2256  (1969). 

98.  Zipf,  E.C.,  Bull.  Am.  Phys,  Soc.  10,  179  (1965). 
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group  has  discussed  these  various  measurements  and  have  suggested  that  incom¬ 
plete  vibrational  relaxation  due  to  low  operating  pressures  (a  few  torr)  may 
explain  the  variety  of  results in  turn,  they  have  derived  their  own 


set  of  pooling  rate  constants  obtained  in  high-pressure  studies  of  e-beam  ex¬ 
cited  Ar/Ng/NO  mixes  (several  thousand  torr).  We  have  chosen  to  use  their 
values  (k12  =  8  x  10"11  cm3/sec  and  k13  ~  2  x  10"11  cm3/sec)^94^ 
as  being  best  representative  of  the  conditions  of  interest  in  this  study. 

lhe  quenching  rate  of  N2(A)  by  N-atoms  is  one  of  the  important  N-(A) 
loss  processes  in  N2  discharges,  lhe  values  of  the  rate  constant  for  this 
process  fall  into  two  groupings:  Young, ^  Wray,^100^  and  Setser^101^ 
have  all  obtained  k17  -  5  x  10"11  cm3/sec  while  Thrush, '102^  et  a!., 

/mM  A/  f 


and  Weinrab, 


(103) 


et  el.  have  obtained  rate  constants  a  factor  of  10  and  100 


times  lower,  respectively.  The  most  definitive  experiment  appears  to  be  that 

of  Setser,  who  directly  monitored  the  time  dependence  of  both  the  N2(A) 

state  loss  and  N(^P)  formation.  These  results  not  only  confirm  the 

Young^")  and  Wray^100^  work  but  also  indicate  that  A  state  quenching  to 

form  the  N(2P)  state  is  the  dominant  loss  channel.  Therefore,  we  have 

-11  3 

chosen  to  use  the  5  x  10  cm  /sec  for  ky  in  these  studies. 

The  N(2P)  states  are  quenched  by  N(  S),  electrons  and  N2(X) 
species.  The  final  products  of  the  quenching  processes  have  not  been  defini¬ 
tively  identified  although  they  must  surely  involve  the  lower  N(20)  and 
N(4S)  states,  we  assume  that  the  ground  state  N(4S)  atoms  is  formed  based 

O 

upon  other  N(  0)  quenching  evidence  which  indicates  that  thp  0-state  quench¬ 


ing  rate  (D>S)  by  N2(X)  is  two  orders  of  magnitude  faster  than  the  P-state 
quenching  rate.  If  P-state  quenching  proceeds  via  the  D-state  intermediate, 
then  the  overall  rate,  P>S,  will  be  limited  by  the  P»D  rate. 


“991  Young,  R.A.  and  St.  John,  G.,  0.  Chem.  Phys.  48,  895  (1968). 

100.  Wray,  K.L.,  J.  Chem.  Phys.  44,  623  (1966). 

101.  Meyer,  J.A.,  Setser,  O.W.  and  Stedman,  D.H.,  J.  Phys.  Chem.  74,  2238 

(1970).  ” 

102.  Thrush,  B.A.,  0.  Chem.  Phys.  47,  3691  (1967). 

103.  Weinreb,  M.  and  Mannella,  G.G.,  J.  Chem.  Phys.  50,  3129  (1969). 

93 


Ihere  are 
cm3/ sec 


10 


-16 


re  four  measurements  of  P-state  quenching  by  N«(X):  k1Q  <  1  x 
,<1M>  <  3  x  10-16  oW1C5)*  6  x  !0-?4  cm3/«cT<106> 

-19  3  xy  T1071 

3  x  10  cm'/sec.  1  Again,  it  is  difficult  to  choose  among  these 
values.  Husain  et  al.^03^  argue  that  Kaufman's  result^06^  is  too  fast 
and  may  reflect  quenching  by  other  transient  species  that  appear  in  their 
higher  power  discharge.  In  addition,  Husain's  result  is  a  direct  time  re¬ 
solved  quenching  measurement  and  should  be  reliable.  We  have  (somewhat  arbi- 
trarily)  chosen  kjo  *  2  x  10"x  cm  /sec  based  upon  the  average  of  the 
two  values  which  are  in  close  agreement. 

P  state  quenching  by  N{  S)  is  important  at  high  energy  loadings  where 
significant  N(4S)  are  produced.  Wayne,  et  al.  reported  a  preliminary  mea¬ 
surement  kx3  *  1  x  10"11  cm3/sec  accurate  to  within  a  factor  of  two.^03^ 
lliis  value  was  used  in  the  present  study. 

An  estimate  of  superelastic  quenching  of  N(  P)  by  electrons  was  ob¬ 
tained  by  using  the  up  excitation  rate  constant  for  (28) 


e$  +  N(4S)  >  N(2P)  +  e$ 


(28) 


and  the  detailed  balance  expression: 

^28  sihx  .-AE/kT 
*20  g(4$) 

where  aE  =  3.58  eV.  Ihe  rate  constant  k,0  has  been  calculated  by  Al i ^ 109) 

co  mol 

who  integrated  the  cross  sections  of  Henry,  et  al.v  ;  for  process  (28) 


104.  Donovan,  ft.J.  and  Gi I lespie, ~H.M. ,  Specialist  Periodical  Reports. 
Reaction  Kinetics  -  Vol  I,  P.G.  Ashmore  Ed.  (the  Chemical  Society 
Burlington  House,  London,  1975)  Chapt.  2. 

105.  Husain,  0.,  Kirsh,  L.J.  and  Wiesenfe’d,  O.R.,  Faraday  Discuss.  Chem. 

Soc.  No.  53,  201  (1972). 

106.  Lin,  C.-L  and  Kaufman,  F. ,  J.  Chem.  Phys.  55,  3760  (1971) ^ 

107.  Noxon,  J.F.,  J.  Chem.  Phys.  36,  926  (1962). 

108.  Golde,  M.F.  and  Ihrush,  8. A.,  Faraday  Discuss.  Chem.  Soc.  No.  53,  233 
(1972). 

109.  Ali,  A.W.,  U.S.  Naval  Res.  Lab  Plasma  Dynamics  lech.  Note  24  (1970);  see 
also  DNA  Handbook,  2nd  Ed.  DNA  1948  H,  p.  20-46. 

110.  Henry,  R.J.W.,  Burke,  P.G.  and  Sinfailam,  A-L.,  Phys.  Rev.  178,  218 

(1969).  — 
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over  a  Maxwell  Boltzmann  electron  energy  distribution.  Using  the  detailed 
balance  expression,  we  obtain  k2Q  ~  1  -  3  x  10  cm  /sec  for  the  elec¬ 
tron  temperature  range  (7fi  ^  1-2  eV)  of  interest  here.  It  should  be  remem¬ 
bered  that  we  showed  in  Section  IV-C-2  that  the  Maxwel 1-Roltzmann  electron 
energy  distribution  is  only  a  crude  approximation  to  the  true  energy  distribu¬ 
tion.  Iherefore,  we  expect  all  of  our  superelastic  quenching  rate  constants 
to  be  somewhat  in  error,  probably  on  the  low  side.  However,  the  effect  of  the 
errur  on  the  kinetics  is  small  because  process  (20)  is  not  the  dominant 

p 

N(£P)  quenching  process  due  to  the  low  electron  density  in  the  afterglow. 

We  also  note  that  the  superelastic  quenching  of  N2(A)  by  electrons  was 

evaluated  in  a  manner  similar  to  that  for  k2Q.  In  this  case,  the  up- 

excitation  rate,  k^,  was  directly  obtained  from  our  Boltzmann  code  results. 

However,  the  use  of  detailed  balancing  still  requires  the  assumption  of  a 

Boltzmann  energy  distribution.  Therefore,  the  caveat  regarding  the  error  in 

-9  3 

k2Q  is  also  applicable  to  our  value  of  k^a,  -  10  cm  / sec. 

We  assume  that  only  N(^S)  atoms  can  recombine  to  form  N2  via  the 
three-body  process  (21),  the  only  net  N(4S)  loss  process  in  our  model.  We 
use  the  Leeds'  group  recommended  value  k21  «  8.3  x  10-J^  exp  (500/ T) 
cm6/sec/111^  Other  recombination  processes  such  as: 

N(  i)  +  fO) - ^  N2*  +  N2  (29a) 

NfS)  +  fO) - N2*  +  N2  (29b) 

are  probably  dominated  by  quenching  processes,  (e.g.,  process  (18))  especially 
considering  that  (29a)  and  (29b)  must  remove  even  more  energy  than  (21).  This 
suggests  the  k29a  and  k2gb  will  be  even  slower  than  k21  and  both  pro¬ 
cesses  can  be  neglected. 

Finally,  we  note  that  the  following  other  processes  can  also  be 
neglected  because  they  are  much  too  slow  to  affect  significantly  any  of  the 
kinetics: 


111. Baulch,  D  l.,  Orysdale,  0.0.  and  Horne,  D.G.,  Evaluated  Kinetic  Data  for 
High  Temperature  Reactions  -  Vol  II,  (Buttersworth  and  Co.,  London, 

19/3;  p.  25. 
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N2(A)  - 

N2(X)  +  hv 

7  .  2 

(30) 

N(2P)  - 

N(4S)  +  hv 

x  ■  12  sec^*3) 

(31) 

N2(A)  +  n2 

— •-  Products  k 

£  IQ"19  cm3/sec^^ 

(32) 

N2(A)  — 

Diffusion  D  ■ 

0.202  cm2/sec(115) 

(33) 

at  1  atm,  300*K 

4.  Results 

The  set  of  coupled  rate  equations  that  describe  the  21  processes  in  our 

model  ware  numerically  integrated  on  an  IBM  computer  using  the  Gear 

algorithm^1*6’ a  sixth  order  predictor-corrector  method  that  also  allows 

for  variable  integration  order  and  variable  step  size.  This  method  is  usually 

applied  to  stiff  sets  of  differential  equations  in  which  the  characteristic 

time  constants  (eigenvalues)  have  widely  differing  values  as  is  the  case  in 

the  ^(A)  kinetics  code.  For  example,  N,,(A)  state  loss  is  controlled  by 

N(4$)  quer*  »i ng  which  has  a  typical  characteristic  time  of  2  usee  (at  0. 

2  -15  2 

*  10  A/cm  and  E/N  =  1  x  10  V-cm  )  while  the  characteristic  time  for 
N(4S)  three-body  recombination  is  about  460  usee  under  the  same  conditions. 

The  computer  program  was  designed  to  optimize  the  N2(A)  state  produc¬ 
tion  relative  to  atom  production  during  the  e-beam  pulse.  This  was  accom¬ 
plished  by  leaving  the  e-beam  and  discharge  field  on  as  long  as  the  ratio, 

A  O 

N^(A)/(N(  S)  +  N  (  P) ) ,  continued  to  increase.  Once  the  maximum  in  the 

ratio  was  reached  the  e-beam  and  discharge  field  were  shut  off.  Typical  pump- 

2 

ing  times  were  of  the  order  of  1  usee  or  less  for  0  ,  =  10  A/cm  to  as 

GO  o 

much  as  several  hundred  microseconds  for  0  ^  =  10  pA/cmS 


112.  Shemansky,  O.E.,  d.  Chem.  Phys.  51,  689  (1969). 

113.  Garstang,  R.H.,  in  Aurora  and  Airglow,  (Pergamon  Press,  Now  York,  N.Y. 
1956)  p.  324. 

114.  DNA  Reaction  Rate  Handbook,  Second  Edition,  DNA  1948H,  March  1972, 
Chapt.  20. 

115.  Zipf,  E.C.,  J.  Chem.  Phys.  38,  2034  (1963). 

116.  Gear,  C.W.,  Comm.  ACM  14,  176  (1971). 

117.  Gear,  C.W.,  Comm.  ACM  14,  185  (1971). 
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In  order  to  check  that  instabilities  in  the  calculation  were  not  being 

generated  by  the  abrupt  e-beam  turnouff,  several  comparisons  were  made  between 

a  continuous  integration  scheme  (from  10  sec  to  10  sec)  and  a  two- 

step  integration  scheme  that  first  integrates  up  to  the  e-beam  shutoff  point 

and  then  uses  the  concentrations  at  the  point  as  initial  conditions  for  a  sec- 

-2 

ond  integration  from  e-oeam  shutoff  to  *0  sec.  In  all  cases  it  was  shown 
that  either  method  gave  the  same  results  within  the  error  limit  of  the  code. 
This  relative  error  was  typically  set  at  10"6  of  the  maximum  value  of  any 
species  concentration. 

Figures  31  through  36  display  the  computer  output  in  graphical  form  for 
different  values  of  Jeb  rr*  E/N.  Of  the  s^ven  species  represented  by  the 
model,  four  were  chosen  to  be  displayed:  N,(A)  (closed  circles),  secondary 
electrons  (open  squares),  NrS)  (open  circles)  and  N(  P)  (closed  squares). 

In  general,  it  is  seen  that  the  temporal  behavior  of  a  given  species  is 
approximately  the  same  irrespective  of  the  value  Jgb  or  E/N.  Electron  den¬ 
sities  rise  at  a  linear  rate  proportional  to  the  source  function  except  at 
high  E/N  and  low  where  secondary  ionization  is  important.  In  this  case, 
exponential  growth  is  observed.  Steady-state  electron  densities,  ne,  are 
reached  if  the  e-beam  is  left  on  for  a  time  long  compared  with  the  recombina¬ 
tion  time.  Under  no  conditions  did  we  ever  observe  a  catastrophic  growth  of 
electron  density  as  a  result  of  a  volumetric  instability.  While  the  electron 
density  grows  linearly  in  cime,  both  N(^S)  and  Ng( A)  grow  quadratically  as 
both  are  produced  by  secondary  electron  impact  on  ground  state  Ng,  i.e., 
dN2(A)  «dN(4S)  <*nedt  tdt.  Likewise,  N(2P)  grows  as  t5  since  it  is 
produced  by  N(^S)  quenching  N2(A),  i.e.,  dN( 2P)  «  [N2(A)][N(4S)]  dt  «  t2t2dt. 

Once  the  e-beam  and  discharge  field  are  shut  off  the  electron  density 
decays  at  the  rate  characteristic  of  electron  ion  recombination.  This  rate 
constant  was  taken  to  be  the  room  temperature  value,  k3  (Tg  =  Tg  =  300*K) 
for  all  the  conditions  studied  here.  In  effect,  we  are  assuming  instantaneous 
electron  cooling  down  to  ambient  temperature  and  negligible  thermal  heating. 
The  time  dependence  of  the  recombination  process  is  given  by 

ne{t)  =  neo/(1  +  k3  neo  tJ 
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TIME  (sec) 


Figure  31  ThecResults  of  the  N?(A)  Kinetics  Code  at  Jeh 
io-5  arops/cmz,  E/N  »  4.1  x  10-1°  v-cmz,  T  *  2 
1  atm  N, 


CONCENTRATIONS  (cm" 


CONCENTRATIONS  (cm*3) 


CONCENTRATIONS  ( cm*3) 


Figure  34  The  Results  of  the  N?  (A)  Kinetics  Code  at  Jeb  « 

10"3  amps/cm*,  E/N  «  1.02  x  10-*3  v-cmS  T  >  298*K  and 
1  atm  N2 
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Figure  36  The  Results  of  the  Ng  (A)  Kinetics  Code  at  Jeh  > 

10_1  amps/cm2,  E/N  ■  1.02  x  10“15  v-cm2,  T  ■  298*K  and 
1  atm  Ng 


'3 


where  nan  is  the  electron  density  when  the  e-beam  is  shut  off.  This  leads 

eo  l  i 

to  an  approximate  t  time  dependence  for  t  »  (k^  neQ)  as  can  be 
seen  in  Figures  31  through  36. 

A  key  feature  of  the  graphical  output  is  the  rapid  rate  of  loss  of 
Ng( A)  states  that  typically  occurs  on  the  one-half  to  several  millisecond 
time  scale.  It  is  apparent  that  appreciable  Ng(A)  state  density  can  only  be 
obtained  if  A  states  are  rapidly  extracted  from  the  discharge  region  on  this 
time  scale  or  faster.  Otherwise,  quenching  will  reduce  A  state  population  to 
a  negligible  value.  In  Section  IV-C-6  we  describe  one  method  for  extracting 
the  A  states. 

The  rate  of  loss  of  Ng(A)  states  is  controlled  primarily  by  the  N(  S) 
quenching  process  (17)  whereas  metastable  ionization  loss  processes  (14),  (15) 
and  (16)  and  the  superelastic  relaxation  process  (17a)  play  only  negligible 
roles  in  the  afterglow.  The  A  state  loss  rate  depends  upon  three  quantities: 
(1)  the  absolute  atom  density  at  the  time  of  e-beam  shutoff,  (2)  the  residence 
time,  tr,  in  the  discharge  zone  and  (3)  the  rate  at  which  N(  P)  is 
quenched  to  H(4S).  Assuming  for  the  moment  one  hypothetical  limiting  case 
where  the  h(2P)  to  N(4S)  quenching  rate  is  infinitely  fast  leads  to  an 
approximate  expression  for  the  time  dependence  of  N2(A)  in  the  afterglow: 

m  m  -kNt„ 

N2H(t)  -  N2m(A)  e  r 


where  hg  (A)  is  the  A-state  number  density  at  shutoff;  k  «  k^;  and  N  is 
the  sum  of  hH(^S)  pulse  f/^P)  at  shutoff.  This  expression  describes 
a  situation  in  which  the  ground  state  atom  density  can  be  recycled  (and  is 
effectively  constant).  The  P  atoms,  formed  in  the  A  state  quenching  pro¬ 
cess,  are  being  instantly  converted  back  to  AS  atoms  which,  in  turn,  produce 
more  P  states  through  A  state  quenching.  On  the  other  hand,  if  we  assume  a 
negligible  N(2P)  to  h(4S)  quenching  rate,  then  the  N(4S)  population  can 


be  rapidly  used  up  and  recycling  of  ground-state  atoms  becomes  impossible, 
this  case,  [N^  (A)  -  fZ^-S)]  net  A  states  will  be  left  at  t  ,  «•  if 
N2  (A)  2  hM(4S).  If  NM(4S)  2  Ng  then  A  states  will  be 
lost. 

In  reality,  the  N(2P)  to  N(AS)  quenching  rate  is  sufficiently  rapid 


that  significant  recycling  of  ground-state  atoms  occurs.  Two  processes  domi¬ 
nate  the  2P  quenching:  relaxation  by  N(AS),  process  (18),  and  relaxation 
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by  molecular  N2»  process  (19).  Of  these  two,  process  (19)  dominates  except 
when  the  N(  S)  density  reaches  ~  5  x  10^  cm”^.  Ihis  occurs  for  the 
cases  where  Jeb  j>  10“*  A/cm^.  Since  the  rate  is  at  least 

5000  sec  *  (  *  200  gsec),  each  jtate  can  be  recycled  at  least  five 
times  assuming  a  1  msec  residence  time. 

As  a  result  of  these  considerations  7able  4  presents  the  maximum  number 
of  A  states,  N2  (A),  produced  in  the  discharge  and  its  lifetime  at  a  variety 
of  Jeb  and  E/N  values  at  298*K  and  1  amagat  N2.  We  have  also  tabulated 
the  discharge  energy  deposited  in  the  gas,  Eq,  and  the  fraction,  fM,  of 
discharge  energy  converted  to  metastable  A  state  excitation.  7he  former  quan¬ 
tity  is  given  by  the  expression: 


E, 


(I) 

e(E/N)N  /  idt  -  e(E/N)  NVdneM  eb/2 


where  N  is  the  N2  density  at  one  amagat,  Vd  is  the  electron  drift  velocity 
calculated  from  the  Boltzmann  code,  n^  is  the  maximum  secondary  electron 
density  at  e-beam  shutoff,  *eb  is  the  pulse  time  of  the  e-beam  and  discharge 
field,  and  the  factor  of  two  accounts  for  the  assumed  triangular  time  depend¬ 
ence  of  the  discharge  current,  i.  7he  quantity,  fM,  is  given  by: 

fM  -  (n2m(a)  ae/ed)  x  100 

where  aE  «  9.87  x  10“^  J.  7he  e-beam  energy  deposited  in  the  gas  is  negli¬ 
gible  compared  to  the  discharge  energy  and  is  neglected  in  calculating  fM. 

7he  graphical  output  data  and  lable  4  indicate  that  the  lifetimes  are 
longer  at  low  E/N  and  low  Jgb  because  the  absolute  N(^S)  density  is  lower 

for  these  conditions.  On  the  other  hand,  it  is  seen  that  the  largest  values 
M 

of  N2  (A)  are  obtained  at  high  E/N  and  Jgb.  Iherefore,  some  trade-offs 
will  be  required  to  obtain  optimal  A  state  density  and  flux.  It  should  be 
noted  that  these  lifetimes  were  obtained  from  the  first  e-fold  times  for 
A  state  loss.  For  the  cases  where  N(^S)  is  growing  rapidly  over  this  time 
scale,  A  state  decay  is  not  described  by  a  single  characteristic  time.  In¬ 
stead,  the  temporal  behavior  indicated  in  Figures  31  through  36  should  be  used. 

One  of  the  key  unknowns  in  the  A  state  kinetics  calculations  is  the 
effects  of  gas  heating  (both  vibrational  and  translational).  7he  kinetics 


105 


KINETICS  CODE  RESULTS 


code  described  above  assumes  that  the  gas  is  always  at  298*K  and  ignores  any 
excited  vibrational  state  effects  on  the  rate  constants.  In  addition,  the 
code  does  not  keep  track  of  the  gas  temperature  rise  nor  are  temperature  de¬ 
pendent  rate  constants  included.  Ihese  effects  may  be  important  for  many  of 
the  cases  shown  in  lable  4  where  the  energy  deposition,  EQ,  will  lead  to  as 
much  as  a  500-1000*K  temperature  rise.  The  specific  heat  of  N2  is  about  1 
J/l-atm-'K.  Ihe  uncertainty  arises  because  neither  the  rate  of  heating  (at 
the  E/N  values  used  here)  nor  the  temperature  dependences  of  the  important 
quenching  processes  are  known.  It  is  reasonable  to  expect  that  the  heating 
pulse  evolves  on  a  time  scale  comparable  to  the  time  scale  for  A  state  extrac¬ 
tion  (~  msec).  However,  until  more  information  becomes  available,  this  uncer¬ 
tainty  will  remain. 

Finally,  we  comment  on  the  apparent  population  inversion  of  N(2P)  with 
respect  to  N(4S)  that  occurs  on  the  post  e-beam  time  scale  (see  Figures  31 
through  36).  Ihis  effect  is  seen  to  occur  whenever  the  e-beam  production  of 
N2(A)  exceeds  that  of  N(4S).  The  inversion  is  created  because  the  loss  of 
N(4S)  via  N2(A)  quenching  is  faster  than  N(2P)  quenching.  The  possibil¬ 
ity  of  a  storage  laser  at  3467  A  should  be  considered. 

5.  Effects  on  N^Gas  Density,  Initial  Gas  Temperature  and  Added  Species 

We  discuss  briefly  possible  modifications  in  the  above  model  that  could 
improve  A  state  production  by  reducing  the  rate  at  which  A  states  are  quenched 
by  ground  state  N  atoms. 

A  substantial  reduction  in  the  initial  N?  gas  density  below  one  amagat 

2^4 

would  help  by  slowing  down  the  rate  of  N(  P)  »  N(  S)  relaxation.  This 
would  effectively  lower  the  number  of  times  that  a  N( 4S)  atom  quencher  could 
be  recycled  given  a  fixed  residence  time  in  the  discharge  zone.  The  limit  on 
N2  density  reduction  would  occur  when  N(2P)  relaxation  becomes  dominated 
by  N  atom  quenching  (N(^S)  >  5  x  10^  cm"2).  One  deleterious  effect  of 

~  4 

lowering  the  N 2  density  would  be  the  reduced  rate  of  N(  S)  three-body  re¬ 
combination.  However,  this  would  be  negligible  at  low  0eb  values  where  N 
atom  loss  through  recombination  is  very  slow  anyway. 

A  substantial  reduction  in  the  initial  N2  gas  temperature  would  accel¬ 
erate  the  rate  of  N  atom  loss  via  the  three-body  recombination  process  (21). 
The  best  value  for  the  temperature  dependence  of  the  N(4S)  recombination 
rate  constant  is  given  by  the  Leeds  value^*^ 
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k?1(7j  »  8.3  x  10~34  exp  (500/7)  cm6/sec 

for  100  K  £  T  £  6C0*K  and  M  »  N2.  Ihis  implies  that  a  thirty-fold  faster 
recombination  time  might  be  realized  by  dropping  the  ambient  gas  temperature 
from  300*K  to  100*K.  Unfortunately,  this  reduction  would  have  the  greatest 
impact  for  the  high  energy  loading  cases  where  the  temperature  rise  from  the 
discharge  heating  processes  might  wipe  out  any  beneficial  effects  that  accrue 
from  the  initial  gas  cooling. 

Ideally,  it  would  be  advantageous  to  add  small  concentrations  of  another 
species  which  wou'io  chemically  remove  the  N  atoms.  This  species  would  have  to 
have  two  properties:  (1)  it  must  react  rapidly  with  N(4S)  and/or  N(3P), 
and  (2)  both  the  reactant  species  and  any  products  ot  the  reaction  must  be  in¬ 
efficient  A  state  quenchers.  The  second  requirement  is  the  most  difficult  to 
achieve  since  it  demands  that  reactants  and  products  have  all  their  electronic 
energy  levels  above  the  A-state  or  else,  if  there  are  levels  at  or  below  this 
level,  they  must  be  at  least  as  metastable  as  N2(A).  At  present  we  have  not 
identified  any  likely  candidates. 

6.  Extraction  of  Metastable  N„(A)  from  the  Discharqe  Zone 

- - *'  l 

N2(A)  formed  in  the  e-beam  sustained  discharge  region  could  be  rapidly 
extracted  from  the  discharge  zcne  and  formed  into  a  supersonic  molecular  beam 
using  the  nozzle  sampling  arrangement  shown  schematically  in  Figure  37.  The 
high  throughput  pump  is  used  to  generate  a  high  flow  speed  in  the  discharge 
zone  and  to  pump  most  of  the  mass  flow.  Ihe  nozzle  then  serves  to  sample  the 
flow  at  a  modest  mass  flow  rate  and  form  the  molecular  beam. 

Our  model  calculations  suggest  that  the  lifetime  of  the  N2(A)  states 
is  of  the  order  of  1  msec  which  means  that  the  flow  time  from  the  discharge  to 
Me  point  A  beyond  the  skimmer  whore  the  metastables  will  be  used  must  also  be 
of  that  order  or  less.  For  L  *  10  cm,  a  flow  speed  of  104  cm/sec  (Mach 
number  a  0.28  at  300*K)  is  required  Shock  wave  formation  within  the  dis¬ 
charge  zone  can  be  reasonably  ignored  at  this  flow  speed.  Assuming  that  the 
discharge  cross-sectional  area  is  about  6  cm3,  (comparable  to  the  dimen¬ 
sions  in  the  AERL  small  scale  "minibang"  e-beam  sustained  discharge  device), 
the  total  mass  flow  rate  is  equal  to  60  std.  1/sec  (127  std.  cfm).  Using  this 
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Figure  37  Schematic  Diagram  of  E-Beam  Sustained  Discharge  Flow  System  for 
N2  (A)  Production 
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value  and  the  equation  of  continuity,  we  can  compute  the  area  A*,  of  the  en¬ 
trance  hole  to  the  high  throughput  pump  from:^*^ 


A*  p 


o 


(y*1)/2(t-1) 


60  std.  1/sec 


where  pQ  is  the  number  density  at  the  stagnation  pressure  (PQ  ■  1  atm)  and 
stagnation  temperature  {7_  ■  298*K),  M  is  the  molecular  weight  of  N?  and 

°  p  c 

v  ■  Cp/Cv  for  Ng  «  1.4.  Evaluating  the  above  yields  A*  ■  2.9  cm  . 

Ihis  ‘•e.'ult  is  in  error  by  a  few  percent  given  that  the  stagnation  pressure 
and  temperature  are  about  5  percent  higher  than  the  values  used  above  when  the 
flow  speed  is  Mach  0.28.  Ihe  area  of  the  sampling  orifice,  AN,  is  then 
arbitrary  and  can  be  sized  to  match  the  appropriate  pump  displacement  that 
evacuates  the  skimmer  zone. 


118.  Hill,  P.6,  and  Peterson,  C.R.,  Mechanics  and  Ihermodynamics  of 
Propulsion,  (Addison-Wesley,  Reading,  Mass.,  1965)  p.  4/. 
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0.  METASTABLE  PRODUCTION  FROM  MAGNETIC  ANNULAR  ARC 

Many  atomic  or  molecular  nmtastable  production  techniques  rely  on  elec¬ 
tron  impact  excitation  processes  to  access  metastable  levels  using  either 
e-beams  or  electrical  discharge  sources,  lhese  methods  usually  require  a  high 
average  electron  energy  in  order  to  achieve  efficient  pumping.  However,  hot 
electrons  can  also  be  detrimental  to  the  metastable  production  efficiency  if 
ionization  of  the  metastable  can  occur.  For  example,  in  helium,  the  23$ 
metastable  level  (~  20  eV  above  the  ground  state)  is  within  5  eV  of  the  ioni¬ 
zation  continuum.  The  cross  section  for  the  ionization  process 

He(23S)  +  e  >  h*e+  +  e  (34) 

is  about  two  orders  of  magnitude  larger  than  the  cross  section  for  the  excita¬ 
tion  process 


He(l  lSQ)  +  e  >  He(23S)  +  e 


(35) 


at  the  electron  energy  which  maximizes  process  (35).  It  is  obvious,  then, 
that  metastable  ionization  can  severely  limit  the  ratio  of  metastable  to 
ground-state  densities.  An  alternative  production  method  is  described  below 
which  utilizes  metastable  formation  from  three-body  atomic  ion-electron  re¬ 
combination  processes,  which  are  most  efficient  at  electron  energies  much 
below  that  required  for  metastable  ionization. 

We  begin  the  discussion  with  a  brief  overview  of  the  salient  features  of 
this  method  and  then  describe  in  more  detail  various  aspects  of  the  ionization 
source,  electron  cooling  processes  and  the  recombination  processes. 

The  basic  idea  is  to  create  a  highly  ionized  plasma  in  which  the  elec¬ 
tron  temperature  is  initially  high  and  the  heavy  particle  temperature  is  low. 
The  plasma  is  allowed  to  expand  into  a  lower  pressure  stagnation  region  where 
the  electrons  cool  via  collisions  with  cold  heavy  ions  and  neutrals.  If  the 
electrons  cool  sufficiently  (down  to  several  tenths  of  an  eV),  volumetric 
three-body  atomic  ion-electron  recombination  of  the  type 
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:  iu  occur  on  a  timescale  fast  compared  to  ion  diffusion  processes  to  the 
walls.  In  process  (36)  He**  represents  any  helium  excited  state.  These 
highly  excited  neutrals  will  cascade  down  via  a  combination  of  radiative  and 
collisiona!  de-excitation  processes  eventually  terminating  in  the  low  lying 
metastable  states. 

The  plasma  source  is  a  magnetic  annular  arc  (MAARC)  also  known  as  a 
plasma  jet  or  magneto-plasma  dynamic  arc  jet.  This  source  is  characterized  by 
a  high  degree  of  ionization  as  calculated  from  the  Saha  equation:^*9^ 


r  9  + 

N  m  [N  A  -~ 
e  1  a  g. 


~mTe 


where 

-3 

N„*  electron  number  density  (cm  ) 

*3 

N  «  ground-state  atom  density  (cm  ) 

A  «  6.04  x  1021  cm**3  -  eV~3/2 

Te«  electron  temperature  (*K) 

I  -  ionization  potential  (ev) 

g+*  degeneracy  of  the  Ion  state 

ga*  degeneracy  of  the  atomic  ground  state 

For  example,  10  percent  ionization  in  helium  (I.P.  »  24.5  eV)  can  be  achieved 

1C  *5 

with  T  -  1-16  eV  for  a  helium  density  of  10  cm  .  However,  it  must 
be  remembered  that  the  Saha  equation  is  applicable  only  for  Boltzmann  equilib¬ 
rium  distributions.  The  applicability  of  (37)  to  the  present  disequilibrium 
situation,  in  which  a  high  electron  temperature  and  low  gas  temperature  are 
required,  demands  that:  (1)  the  time  scale  for  electron-electron 
equilibration,  ee,  be  short  compared  to  the  residence  time  in  the  arc 
region,  and  (2)  the  electron-ion  equilibration  time  scale,  be  long 
compared  to  the  resi-  dence  time  in  the  arc,  i.e., 


(36a) 


T 


ee 


1 


NQ 
eHee 


7T 


<  tFLOW 


T 


el 


N^e 


mc  >  T  FLOW 


(38b) 


where 

p 

Qee  -  the  cross  section  for  electron-electron  scattering  (cm) 

Qd  ■  the  cross  section  for  electron-ion  scattering  (cm2) 

ce  ■  average  electron  speed;  1/2  mec2  ■  3/2  kTe 

Eq.  (38a)  insures  that  there  is  a  sufficient  time  to  build  up  the  high  energy 
tail  of  the  Boltzmann  distribution  which  is  the  source  of  the  ionization  and, 
(38b)  insures  that  ion  heating  by  the  electrons  is  slow,  lhe  cross 
sections,  Qfie  and  QeJ  can  be  obtained  from  the  coulomb  scattering  formula: ^12(^ 


4.4  x  10' 


-14  ,4 


T„ 


l"Al 


4.4  x  10 


-13 


cm 


(39) 
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for  singly  charged  ions  (Z  ■  1).  Assuming  N  «  10  cm  and  T  *  1  eV 

-9  -5e  e5 

yields  Tee  *  4  x  10  sec  and  =  3  x  10  sec.  Given  the  10  cm/sec 

flow  speed  in  the  arc  region  implies  that  a  flow  dimension  Lj  <  3  cm  would 

be  required  to  satisfy  (38a)  and  (38b). 

A  schematic  diagram  of  the  source  is  displayed  in  Figure  38.  The  device 

consists  of  a  coaxial  pair  of  electrodes,  the  cathode  in  the  center  and  an 

anode  which  can  also  serve  as  a  nozzle.  The  electrons  are  heated  in  the  arc 

region  over  a  distance  L^,  and  then  are  cooled  as  they  traverse  the  distance 

l£.  The  B-field  is  used  to  stabilize  the  discharge  and  provide  a  means  for 


120.  Spitzer,  L.,  Physics  of  Fully  Ionized  Gases,  2nd  Edition,  Interscience 
Publishers  (Wiley),  New  York,  1962. 
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achieving  a  relatively  large  voltage  in  the  plasma  stream  even  when  the  plasma 
has  a  large  conductivity,  7he  criteria  for  achieving  this  impedance  is  to 
make  wexe  >  1  where  «e  -  eB/cme  is  the  gyro  frequency  of  the  electrons 
and  te  ■  l/NeQeIce  is  the  electron  ion  collision  time  for  electrons 
which  is  appropriate  for  highly  ionized  plasmas.  Note  that  xe  is  much 
smaller  than  the  ion  heating  time,  xgj.  Since  the  magnetic  field  is  aligned 
essentially  parallel  to  the  arc  channel  walls,  and  since  the  electrons  undergo 
complete  gyro  orbits  around  the  field  lines  between  collisions,  their  mobil¬ 
ity  perpendicular  to  the  field  is  reduced.  Ihe  required  B-field  is  then  given 
by: 


cm. 


8  >  if  NeQeIce 


(40) 
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For  1  eV  average  energy  electrons  in  helium,  B  >  14  G  at  Ng  *  10  cm  . 

Ihe  MAARC  operating  voltage  can  be  expressed  by  the  following  empirical 
formula:^121) 


V  •  V0  +  UCB* 


(41) 


where 


0  2 
1/2  Mjlr 

B  « 

£  *: 


Electrode  voltage  drop  -  30  V  for  helium  =  constant. 

=  E.  .  . .  -  25  eV 

ionization  - 

Magnetic  field  strength  (Gauss) 
radial  separation  between  electrodes  (cm) 


2 

Ihe  quantity,  1/2  M.lr,  is  the  ion  kinetic  energy  set  equal  to  the  ion- 

*  C  e 

ization  energy.  For  helium,  the  value  of  Uc  is  equal  to  ~  5  x  10  cm/sec. 
Hence,  the  U  Bfc  voltage  drop  ~  7.5  V  when  B  *  300  G  and  i  *  1/2  cm.  The  in- 

V  ' 

put  power  to  the  plasma  is  given  by: 


121. Patrick,  R.M.  and  Schneiderman,  A.M.,  Jour.  AIAA  £,  283  (1966). 
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piN  "  ^o^particle 

-  (1015  cm"3)  (10s  cm/sec)  (1.6  x  10"19  (1  eV)  -  16  W 


The  MAARC  operating  current,  I,  is  given  by 

I  ■  16  W/7.b  V  ■  2.1  amps 

The  electrons  will  cool  as  they  expand  out  of  the  arc  (Region  I)  into 
the  recombination  zone  (Region  II).  Ihe  main  cooling  processes  are  momentum 
transfer  collisions  to  the  cold  ions,  and,  subsequently,  the  ions  will  loose 
their  energy  by  symmetrical  charge  transfer  reactions  to  the  background  neu¬ 
trals.  For  efficient  cooling,  the  time  scales  for  electron-ion  momentum 
transfer,  ej,  and  charge  exchange  cx  must  be  less  than  the  flow  time. 

We  have  previously  shown  that  2  3  x  sec  and 

^cx  1 

TCX  “Cf  “ 

»  10  psec 

Therefore,  reI  and  tcx  <  TpLQW  can  be  satisfied  if  the  flow  time  is  of 

the  order  of  100  nsec.  This  requires  that  for  Region  II  L„  ~  10  cm  in  the 

5  ^ 

flow  dimension  for  a  flow  speed  of  10  cm/sec. 

Table  5  indicates  the  extent  of  "ideal"  electron  cooling  assuming  equi- 
partition  of  the  initial  electron  energy  as  calculated  for  several  ionization 
levels.  This  is  obviously  only  a  lower  limit  to  the  true  electron  temperature 
in  the  recombination  zone  since  the  electrons  can  be  heated  up  again  using  the 
energy  liberated  from  the  recombination  processes  which  have  been  neglected  in 
the  calculations  displayed  here.  The  first  column  in  Table  5  indicates  the 
level  of  fraction  ionization,  the  second  column  indicates  the  electron  tem¬ 
perature,  7.  in  the  arc  required  by  the  Saha  equation  to  reach  this  ioniza- 

“  ■  ICO 

tion  level  at  the  neutral  density  N  =  10  cm  ;  the  third  column  is 

a 

the  final  electron  temperature  in  the  recombination  zone  obtained  from 
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TABLE  5 


EQU1PARTITI0N  CALCULATION 


Nei/Na 

Tei(eV) 

V(eV) 

o(cm®/sec) 

Nei 

'  ( S2C) 

1.0 

1.45 

0.725 

2.3  x  lO-H 

1  x  1015 

43  x  10"6 

0.5 

1.35 

0.47 

8.4  x  lO"11 

5  x  loM 

24  x  10-6 

0.1 

1.16 

0.13 

5.4  x  lO"9 

1  x  1014 

2  x  lO"6 

0.01 

0.96 

~  0.035 

2  x  10*7 

1  x  1013 

0.5  x  lO"6 

117 


(43) 


2  Nek7ei  +  2  Nak^a  "2  (Ne  +  Na^^ef 


where  kTa  is  the  room  temperature  thermal  energy  of  the  ions  and  neutrals; 
the  fourth  column  indicates  the  overall  recombination  rate  constant  at  tem¬ 
perature  using  the  Hinnov-Hirschberg  formula:^^ 


a 


-27 

5.6  X  10  4 


ana  the  sixth  column  is  the  characteristic  time  for  recombination  at  T  ^  and 
the  initial  electron  density  Nei-  (column  5).  From  the  data  in  Table  5  it 
can  be  seen  that  the  electrons  will  cool  sufficiently  such  that  the  re¬ 
combination  time  is  always  commensurate  with  or  faster  than  the  flow  time 
(-  100  nsec).  Iherefore,  it  is  possible  to  generate  as  much  as  100  percent 
ior.i2ation  in  the  arc  and  still  cool  the  electrons  sufficiently  such  that 
efficient  volumetric  recombination  occurs. 

Ihe  time  evolution  of  the  metastable  neutral  states  generated  from  the 
recombination  of  atomic  ions  and  electrons  is  a  complicated  dynamical  process 
which  depends  critically  on  the  competition  between  a  number  of  processes. 
Although  it  is  reasonable  to  expect  that  the  ratio  of  singlet  to  triplet  ex¬ 
cited  states  will  be  -  1/3,  the  absolute  fraction  of  metastable  atoms  cannot 
be  estimated  without  extensive  modeling.  7o  amplify  this  point  note  that 
Table  5  indicates  that  the  recombination  time  at  high  fractional  ionization 
levels  is  commensurate  with  the  electron  cooling  rates  discussed  previously. 
In  addition  both  superelastic  processes  of  the  type 

He(23S)  +  e  >  He^S)  +  e  (44) 

and  spin  exchange  inelastic  processes  of  the  type 

He(23S)  +  e  »  He(21P)  +  e  (45) 


122.  Hinnov,  E.  and  Hirschberg,  J.G.,  Phys.  Rev.  125,  795  (1962). 
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peak  at  low  energies  and  thus  compete  with  the  desired  recombination  process. 
Process  (44)  destroys  metastables  directly  and  process  (45)  leads  to  meta¬ 
stable  destruction  via  the  optically  allowed  transition  He(2*P  *  1*S). 

Rather  than  perform  the  necessary  modeling  calculations,  we  have  at¬ 
tempted  to  gain  insight  into  metastable  production  yields  using  the  experi¬ 
mental  data  of  the  Princeton  stellarator  group.  These  experiments  havn 

investigated  electron  loss  rates,  electron  cooling  rates,  electron  and  excited 
state  densities  in  the  afterglow  of  a  highly  ionized  plasma  produced  by  the 
pulsed  stellarator  source.  The  range  of  electron  densities 
(1011  <  Ne  <  3  x  1013  cm”3)  and  electron  temperatures 
(0.04  eV  <  Te  <  1  eV)  correspond  to  the  conditions  we  expect  in  the  MAARC 
generator.  Although  the  main  thrust  of  the  Princeton  work  was  involved  with 
time  resolved  population  densities  in  highly  excited  states  (principal  quantum 
number  >  5),  some  measurements  for  the  n  «  2,  3  and  4  levels  have  been 
obtained.  These  measurements  indicate  that  although  the  highly  excited  states 
are  in  Boltzmann  equilibrium  at  the  electron  temperature,  the  populations  of 
the  low  lying  states  (n  ■  2,3  and  4)  are  orders  of  magnitude  below  what  would 
be  expected  based  upon  the  Boltzmann  factor  for  the  high  states.  Using 
published 

results^3^  we  estimate  a  fractional  metastable  production  -  10“6  -  10“^  of 
the  ground-state  density  for  initial  helium  densities  1. 5-3.0  x  1013  cm-3  and 
electron  temperatures  ~  0.2-0. 4  eV. 

In  conclusion,  it  appears  as  though  the  metastable  loss  processes  such 
as  (44)  and  (45)  are  much  faster  than  the  processes  which  populate  the  meta¬ 
stable  levels  from  the  higher  lying  states,  and  that  the  dominant  product  from 
recombining  atomic  plasmas  is  the  ground  state  of  the  atom. 


123.  Johnson,  L.C.  and  Hinnov,  E.,  Phys.  Rev.  _187,  143  (19G9). 
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V.  STORAGE  OF  KETAS TABLES 


One  can  store  metastables  in  a  volume  confined  by  material  surfaces  for 
a  length  of  time  which  varies  inversely  with  (a)  the  rate  at  which  the  meta¬ 
stables  collide  with  the  surfaces,  and  (b)  the  probability  of  deactivation  on 
the  surfaces.  Since  the  latter  is  expected  to  decrease  with  decreasing  tem¬ 
perature,  one  is  led  to  consider  storage  in  a  cryogenic  container. 

There  is,  however,  another  more  speculative  approach  to  this  problem. 
Recently,  several  proposals  for  confining  atoms  and  molecules  with  laser  light 
have  appeared  in  the  literature, C 124-126)  Asf,j^n,  for  example,  estimates 
that  106  atoms/sec  can  be  confined  in  a  laser-beam  trap  of  capacity 
-  3  x  10^  atoms;  atoms  whose  velocities  range  between  0  and  3  x  10^  cm/sec 
are  captured  in  the  trap.^126^  Ashkin’s  trap  consists  of  two  counter- 
propagating  TEM^q-  mode  laser  beams  focused  to  spot  sizes  of  several  um 
(Figure  39).  The  focal  spots  are  -  2  cm  apart.  As  shown  in  Figure  <10,  a  re¬ 
quirement  is  that  the  laser  light  be  nearly  resonant  with  an  allowed  transi¬ 
tion  of  the  atoms  which  are  to  be  trapped. 

Two  mechanisms  work  together  to  confine  atoms  in  the  counter-propagating 
beam  configuration  shown  in  Figure  39.  First,  there  is  a  net  longitudinal 

force  on  an  individual  atom  which  tends  to  slow  it  down  if  u  <  uQ.  This 
force  may  be  understood  from  Figure  41.  An  atom  is  continuously  absorbing 
laser  light  and  emitting  it  spontaneously  and  isotropically.  In  the  atomic 
rest  frame,  the  light  wave  opposed  to  the  atomic  motion  appears  shifted  up  in 
frequency,  while  the  other  wave  is  shifted  down  in  frequency.  Since  the 
opposing  beam  is  thus  shifted  closer  to  the  atomic  resonance  frequency  uQ, 
an  atom  moving  to  the  right  absorbs  more  photons  from  it  than  from  the  co¬ 
propagating  light  wave  and  is  slowed  down  by  the  net  transfer  of  momentum. 
Similarly,  an  atom  moving  to  the  left  would  absorb  more  strongly  from  the 
right-going  light  wave,  hence,  it  too  would  be  slowed  down. 


124.  Letokhov,  V.S.,  Minogin,  V.G.  and  Pavlik,  0.,  Opt.  Comm.  19,  72  (1976). 

125.  Letokhov,  V.S.  and  Minogin,  V.G.,  0.  Opt.  Soc.  Am.  69,  413  (1978). 

126.  Ashkin,  A.,  Phys.  Rev.  Lett.  40,  729  (1978). 
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Figure  39  Geometry  of  Ashkin's  Optical  Trap  (Ref.  126).  The  atomic  beam  is 
injected  into  the  interaction  region  through  a  hole  in  the 
mirror  so  that  the  atoms  move  along  the  propagation  direction  of 
the  light.  The  focal  point  of  the* laser  beam  propagating  to  the 
left  (right)  is  labeled  Qj  (Q2).  w  is  the  minimum  spot 
size  of  the  laser  beams.  '  0 
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Figure  40  Energy  Levels  of  an  Atom  Confined  to  an  Optical  Trap.  The 
energy  of  the  upper  level  u  is  T»  0  higher  than  that  of  the 
lower  level  .  The  upper  level  decays  at  the  rate  y-  The 
incident  light  has  a  frequency  u  which  is  smaller  than  u0  fay 
-  50  Y. 
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(a)  LAB  FRAME 

(b)  ATOMIC  REST 
FRAME 

K374« 


V  O) 


w-Aw 


u>  +  Acj 

UO)i  -A/VW 


Figure  41  (a)  An  Atom  moves  to  the  Right  at  Velocity  v  along  the  Direction 

of  Propagation  of  One  of  a  Pair  of  Counter-Propagating  Light 
Beams  of  Frequency 

(b)  Viewed  in  the  Atomic  Rest  Frame,  the  Right-Going  Beam 
Appears  Downshifted,  and  the  Left-Going  Upshifted  in  Frequency. 


Ihe  second  mechanism,  the  ponderomotive  force,  confines  the  atom  in  the 
transverse  direction,  and  may  be  understood  as  follows:  the  root-mean-squars 
(rms)  electric  field  of  the  light  can  exert  an  rms  force  on  the  atom  in  the 
same  way  as  a  static  field.  If  the  rms  field  is  spatially  uniform  this  force 
vanishes,  but  any  nonuniform  field  polarizes  an  atom  and  therefore  exerts  a 
net  "ponderomotive"  force  on  the  atom.  It  may  be  shown  that  for  *  <  the 
sense  of  the  force  is  the  same  as  that  due  to  a  nonuniform  static  field,  the 
atom  being  attracted  to  a  region  of  higher  intensity.  Since  the  laser  in¬ 
tensity  is  highest  along  the  beam  axis,  the  atoms  are  continuously  attracted 
to  this  axis  and  the  ponderomotive  force  exerts  a  stabilizing  effect  on  the 
atomic  confinement. 

We  believe  that  this  technique  may  be  generalized  to  the  storage  of  cer¬ 
tain  metastable-state  atoms  and  molecules.  The  basic  ides  behind  such  storage 
is  that  the  laser  should  be  tuned  to  a  frequency  u  slightly  below  the  fre¬ 
quency  uQ  of  an  allowed  transition  between  the  metastable  state  m  and  an 
excited  state  e  (Figure  42a).  MetastabU-state  atoms  injected  into  the 
optical  trap  are  confined,  while  those  atoms  in  other  states,  such  as  the 
ground  state  g,  pass  through  the  trap  without  interacting  appreciably  with  the 
light.  The  trap  thus  acts  as  an  atomic  filter  which  confines  only  those  atoms 
which  are  in  the  metastable  state.  A  possible  candidate  species  for  such  con- 

3 

finement  is  atomic  He  in  the  2  S  state  (Figure  42b),  which  would  require  a 

i  i 

laser  nearly  resonant  with  either  the  1.08-wm  2  S-2  P  or  the  0.389-pm 

23$-33P  transition.  From  the  work  of  Ashkin/*^  we  estimate  that  the 

2 

required  laser  intensity  would  be  ~  1  W/cm  and  the  required  minimum  laser 
beam  spot  size  ~  10  urn.  A  beam  of  23S  He,  for  example,  might  require  a 
quasi-CW  laser  nearly  resonant  with  either  the  l.OS-um  23S-23P  or  the 
0.389-pm  23S-33P  transition. 
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Figure  42  (a)  Energy  Levels  of  an  Atom  of  Interest,  Showing  the  Ground 

State  g,  Metastable  State  m,  and  Excited  State  e. 

(b)  Energy  Levels  of  He  Corresponding  to  Those  of  (a). 


125 


